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ABSTRACT 

The  utility  of  simple  particle  suspensions  in  modelling  complex 
fluids  has  motivated  considerable  study  in  fundamental  suspension 
behaviour.  Unfortunately,  most  of  these  studies  have  been  restricted 
to  examining  suspension  behaviour  in  simple  shear  flows,  neglecting 
the  behaviour  in  extensional  flows. 

The  orifice  jet  thrust  technique  was  used  to  experimentally 
investigate  the  behaviour  of  fibre  suspensions  in  extensional  flows. 
Fibres  with  axis  ratios  of  85  to  340  were  suspended  in  60%  and  66-2/3% 
sugar  solutions  at  concentration  levels  of  0.1%  and  0.2%.  Thrust 
reduction  was  observed  for  suspensions  with  shear  viscosities  similar 
to  the  shear  viscosity  of  the  suspending  medium.  The  thrust  reduction, 
which  may  be  viewed  as  an  increase  in  elongational  viscosity,  became 
more  pronounced  with  increasing  axis  ratio  and/or  fibre  concentration 
as  predicted  by  the  specific  elongational  viscosity  theory  of 
Takserman-Krozer  and  Ziabicki* 

An  unusual  behaviour,  indicated  by  apparently  steady-state  pres¬ 
sure  and  thrust  readings  suddenly  changing  to  new  steady-state  readings, 
was  observed  during  some  of  the  experiments.  This  occurred  for  both 
the  suspensions  and  the  Newtonian  suspending  medium. 
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CHAPTER  I 
INTRODUCTION 

1.1  GENERAL 

The  rheological  behaviour  of  solid-liquid  suspensions  has 
enjoyed  considerable  engineering  interest  in  a  wide  variety  of 
applications.  One  field  of  major  importance  involves  pipeline  trans¬ 
port  of  solids.  Fairly  recent  endeavours  involving  suspensions  are 
concerned  with  modelling  of  viscoelastic  fluids  (33,34,39,40)  and  of 
blood  (15).  An  important  property  of  many  suspensions  is  the 
phenomenon  of  drag  reduction  in  pipe  flows.  This  phenomenon  alone 
has  promoted  considerable  investigation.  An  extensive  review  of 
this  subject  is  presented  by  Kerekes  (21)  who  also  discusses  the 
probable  mechanisms  for  drag  reduction  in  suspension  flows. 

The  prime  difficulty  in  the  handling  of  suspensions  has  been 
in  uniformly  dispersing  the  solid  particles  (1,21).  Since  common 
industrial  slurries  involve  some  heavy  particulate  (such  as  sand) 
and  a  convenient  suspending  medium  (such  as  water),  sedimentation 
is  an  enormous  problem.  Even  in  cases  where  sedimentation  is  not  a 
serious  problem,  as,  for  example,  in  wood  pulp  flows,  plugging  at 
constrictions  is  common.  From  an  experimental  point  of  view,  since 
most  rheological  studies  involve  bulk  measurements,  particle 
coagulation  (or  concentration  gradients)  nullify  the  often  necessary 
assumption  of  material  homogeniety.  In  this  context,  it  is  not 
surprising  that  many  rheological  studies  of  suspensions  have  not  been 
entirely  successful. 
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Although  drag  reduction  studies  have  been  successful  in  terms 
of  relating  drag  reduction  to  suspension  properties,  the  mechanism  of 
drag  reduction  has  been  a  controversial  subject  (21).  Basically,  the 
controversy  has  been  one  of  evaluating  the  importance  of  the  elastic 
contribution  to  drag  reduction.  Resolution  of  this  problem  is  com¬ 
plicated  by  a  lack  of  sufficient  knowledge  on  the  elastic  properties 
of  suspensions.  Clearly,  a  better  understanding  of  this  problem,  and, 
in  fact,  many  of  the  problems  associated  with  drag  reducing  polymers 
(33,34,39,40)  can,  perhaps,  be  explained  on  the  basis  of  the  elastic 
behaviour  of  suspensions  of  relatively  simple  shapes  (29,33,34,39,40). 

A  particular  limitation  of  most  rheological  measurements  is 
that  extensional  flows  or  "stretching"  deformations  of  the  material 
elements  have  largely  been  ignored  (25).  Rapid  extensional  deforma¬ 
tions,  which  are  often  accompanied  by  shearing  deformations,  are 
widespread  in  practice;  as,  for  example,  in  extrusion,  rolling  or 
pressing  operations,  and  flows  through  contractions  or  around  sub¬ 
merged  objects  (24).  Metzner  (24)  points  out  that  the  extensional 
deformations  in  Newtonian  flows  contribute  relatively  minor  resis¬ 
tance  when  compared  to  the  resistances  proffered  by  the  shearing 
deformations.  Viscoelastic  materials,  on  the  other  hand,  exhibit 
very  large  resistances  to  rapid  stretching  deformations  (23,24, 
25,26,37).  Hence,  it  has  been  suggested  explicitly  by  Metzner  (24) 
and  less  explicitly  by  others  (23,26,37)  that  rapid  extensional  flows 
of  viscoelastic  media  may  be  approximated  by  considering  only  the 
stretching  deformations  and  completely  ignoring  the  shearing 
deformations.  Significantly,  this  results  in  a  tremendous 


-  . 


3 


simplification  in  a  number  of  important  problems  (24) .  In  cases  where 
slow  extensional  deformations  prevail  for  viscoelastics,  or  rapid 
deformations  occur  in  fluids  of  a  low  level  of  elasticity,  it  is 
expected  that  both  stretching  and  shearing  deformations  must  be  con¬ 
sidered  in  a  complete  analysis  of  the  problem. 

Although  viscoelastic  behaviour  has  been  observed  in  shearing 
flows  of  fibrous  suspensions  (7,8,10,29),  the  elastic  effect  has 
generally  been  small  (10)  and  complicated  by  time-dependence  and/or 
network  formation  as  well  as  particle  migration  away  from  a  wall 
(10,29) o  Suspension  behaviour  in  extensional  flows  has  not  yet  been 
experimentally  considered .  The  purpose  of  this  study,  therefore,  is 
to  experimentally  consider  thei  behaviour  of  fibre  suspensions  in 
extensional  flows.  This  analysis  will  benefit  from  a  brief  discus¬ 
sion  on  the  available  rheological  properties  of  homogeneous  suspen¬ 
sions  . 

1.2  RHEOLOGICAL  PROPERTIES  OF  ELONGATED  PARTICLE  SUSPENSIONS 

lc2„l  Viscous  Behaviour 

Most  rheological  studies  of  suspensions  have  dealt  specifically 
with  determining  shear  viscosities .  Predictive  viscosity  theories  for 
suspensions  of  simple  particle  shapes  have  been  given  by,  Einstein 
(11)  for  spheres,  Jeffery  (19)  and  Guth  (29)  for  rigid  ellipsoids  and 
Burgers  (29)  for  rigid  ellipsoids  and  rigid  cylinders.  As  noted  by 
Blakney  (6),  all  of  these  theories  can  be  expressed  by: 


.  u 


1  +  cx0C 


(1.2.1) 
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where  yf  is  the  relative  viscosity,  u  the  viscosity  of  the  suspension, 
p0the  viscosity  of  the  suspending  medium,  C  the  volumetric  concentra¬ 
tion  of  particles  and  aQ  a  function  of  particle  orientation.  Agreement 
with  experiment  has  been  good  for  all  particle  shapes  so  long  as  they 
are  non-deformable  (29)  and  very  dilute  (6), 

Experimentally,  it  has  been  found  that  the  linear  viscosity 
function  given  by  Equation  (1,2,1)  no  longer  applies  when  a  critical 
concentration  limit  has  been  reached.  This  critical  limit  has  been 
defined  by  Mason  (10)  such  that  for 

C(AR)2  <  1.5  (1.2.2) 

where  AR  is  the  axis  ratio  of  the  particle  (ratio  of  length  to 
diameter),  particle  interactions  are  negligible.  The  limit  given  in 
Equation  (1,2,2)  is  defined  by  the  minimum  swept  out  volume  required 
by  a  rotating  elongated  particle  (21),  Beyond  this  limit,  the 
dependence  of  viscosity  on  concentration  and  particle  properties  is 
quite  complex.  However,  three  broad  categories  based  on  fibre 
aggregation  have  been  proposed  (10,21): 

i)  dilute  region  where  single  particles  are  free  to  rotate  at 
will  and  any  collisions  are  random, 

ii)  critical  concentration  region  where  particle-particle 

collisions  are  numerous,  but  network  formation  is  almost 
entirely  absent, 

iii)  concentrated  region  where  particle-particle  contact  is  so 
prevalent  that  a  more  or  less  continuous  network  is  formed. 
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In  the  concentrated  region,  the  suspension  exhibits  distinct  non- 
Newtonian  properties. 

Viscosity  measurements  for  rods  in  the  critical  concentration 
region  have  been  substantial  (6,7,8,21,29),  but  attempts  to  correlate 
these  measurements  by  a  2-parameter  equation  have  had  limited 
success  (29),  The  2-parameter  viscosity  equation  is  given  by 


whe  re 


a  +  an  C 
o 


O' 


(1 .2*3) 


(1  -2.4) 


a1  is  defined  as  the  particle  interaction  coefficient  and  is  the 
specific  viscosity,  Nawab  and  Mason  (29)  found  that  Equation  (1,2,3) 
could  be  successfully  used  to  correlate  experimental  data  so  long  as 
the  particles  remained  rigid.  Their  study  is  of  particular  import¬ 
ance  since  they  found  that  fibre  suspensions  began  to  exhibit  a  shear¬ 
thinning  behaviour  when  fibre  deformations  became  apparent,  A 
qualitative  summary  of  the  observations  of  Nawab  and  Mason  (29)  is 
given  in  Table  (1„1).  It  should  be  noted  that  the  available 
viscosity  data  for  suspensions  are,  generally,  from  a  rheological 
standpoint,  limited  due  to  the  narrow  range  of  shear  rates  studied. 

Concentrated  suspensions  have  escaped  reliable  viscosity 
determination  due  to  their  "unusual"  behaviour  of  time  dependence 
and  particle  migration  (10,21,29).  As  a  result,  the  most  successful 
measurements  involving  these  suspensions  have  been  non-rheological  ; 
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for  example,  head-flow  studies,  Bugliarello  and  Daily  (8)  discuss, 
in  detail,  the  behaviour  of  "network"  suspensions  and  conclude  that 
present  continuum  models  do  not  adequately  describe  concentrated 
suspension  behaviour, 

1,2,2  Elastic  Behaviour 

The  two-phase  approach  to  molecular  theories  on  viscoelastic 
behaviour  are  based  on  the  "deformabil ity"  of  a  macromolecul ar 
system.  Elastic  properties  in  the  whole  fluid  can  be  related  to  the 
deformation  and  subsequent  restoration  forces  of  a  macromolecul ar 
system.  The  macromolecule  has  been  considered  from  basically  two 
points  of  view: 

i )  Simple  particle  approach 

The  fluid  is  considered  as  a  suspension  of  simple-shaped, 
elastically  deformable  particles.  When  the  fluid  is  subject  to 
a  stress,  the  particles  deform  and  orient  themselves  in  such  a 
way  that  the  elastic  forces  in  the  particles  and  the  forces 
originating  from  the  orientation  of  deforming  particles  in  the 
suspending  medium,  will  come  to  equilibrium  with  the  external 
stress.  The  bulk  elastic  properties  of  the  whole  fluid  can  be 
described  in  terms  of  the  restoration  forces  within  the 
particles  and  their  subsequent  re-orientation  (12,13,33), 
i i )  Complex  particle  approach 

The  fluid  is  considered  as  a  suspension  of  complex,  chain¬ 
like  particles.  When  the  fluid  is  at  rest,  the  orientation  of 
the  particles,  though  physically  unknown,  can  be  related  to  an 
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equilibrium  of  the  intermolecul ar  and  intramolecular  forces  and 
the  random  forces  of  Brownian  motion.  This  relation  gives  an 
"average"  orientation.  When  the  fluid  is  subject  to  a  stress, 
the  change  in  the  "average"  orientation  can  be  predicted  on  the 
basis  of  the  nature  of  the  stress.  The  bulk  elastic  properties 
of  the  whole  fluid  can  be  described  in  terms  of  the  change  in 
orientation  of  the  macromolecul ar  chain  (27,39,40)  The 
restoration  forces  are  considered  to  be  those  of  brownian 
motion  (18), 

Although  the  relative  merits  of  the  two  approaches  outlined  above 
are  beyond  the  scope  of  this  study,  the  principles  involved  may  be  use¬ 
ful  in  providing  an  explanation  for  the  origin  of  elasticity.  Both 
simple  particle  and  complex  particle  approaches  consider  an  increase 
in  resistance  to  a  stress  due  to  the  presence  of  particles.  This 
increase  in  resistance  can  be  regarded  effectively  as  an  increase  in 
fluid  viscosity  and  is,  by  necessity,  present  in  any  suspension, 
elastic  or  not.  The  elastic  properties  ot  the  fluid,  on  the  other  hand, 
arise  from  the  restoration  forces  which  are  a  consequence  of  an  elastic 
deformation.  At  this  point,  it  is  not  important  whether  the  restora¬ 
tion  forces  are  pre-defined  (simple  particle  approach)  or  related  to 
Brownian  motion  (complex  particle  approach)  only  that  they  are  not 
negiigible  when  compared  to  the  viscous  forces.  It  is  apparent  that 
the  fluid  will  exhibit  only  a  viscous  response  to  a  deformation  if  the 
restoration  forces  are  comparatively  small 

Application  of  the  simple  particle  approach  to  a  dilute  suspen¬ 
sion  of  fibres  appears  to  be  reasonable  enough  since  fibres  are 
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relatively  simple  particles.  Nawab  and  Mason  (29)  indicated  qualitive 
experimental  support  to  this  approach  by  observing  the  "Weissenberg 
Effect"  in  Couette  flows  of  deformable  fibres  (see  Table  1.1). 

Although  the  "WeissenbeYg  Effect"  did  not  occur  until  fibre  deforma¬ 
tions  became  apparent,  the  particles  were  concentrated  enough  to  form 
networks  so  that  it  was  not  possible  to  attribute  the  elasticity 
solely  to  the  recoverable  deformation  of  single  particles.  Indeed, 
the  elastic  counter-rotation  observed  by  Bugliarello  and  Daily  (8) 
was  directly  attributed  to  a  network  relaxation  effect .  Also,  the 
available  data  on  normal  stresses  developed  in  shear  flows  of  fairly 
concentrated  suspensions  (10)  have  shown  the  elastic  effects  to  be 
small .  These  observations  are  consistent  with  the  complex  particle 
approach  in  suggesting  low  levels  of  elasticity  due  to  the  minor 

i 

importance  of  Brownian  motion  on  large  particles.  In  summary,  the 
elastic  effects  exhibited  by  fibre  suspensions  in  simple  shear  flows 
are  generally  small  and  may  be  of  little  importance  in  practical 
situations.  The  larger  effects  expected  in  stretching  deformations 
clearly  underline  the  importance  of  extensional  flows. 

1.2.3  Behaviour  of  Elongated  Particle  Suspensions  in  Extensional  Flows 
Recently,  Takserman-Krozer  and  Ziabkki  (39,40)  analyzed, 
theoretical ly ,  the  behaviour  of  rigid  ellipsoids  in  a  parallel  velocity 
gradient  (extensional  flow).  This  analysis,  which  may  represent  a 
viscoelastic  model  in  steady  extensional  flows,  led  to  the  definition 
of  a  specific  elongational  viscosity 
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(1.2.5) 


(1.2.6) 

n  is  the  elongational  viscosity  of  the  suspension,  nQ  the  elongational 
viscosity  of  the  suspending  medium,  C  the  volumetric  concentration  of 
ellipsoids,  e  the  elongational  deformation  rate,  v the  volume  of  one 
ellipsoid  and  (^—)  the  average  rate  of  energy  dissipation  per  particle, 
which  is  determined  from 

($f)  »  f(AR,a,  nQ)  (1.2.7) 

where 


and  D  is  the  rotational  diffusion  constant  (34).  The  intrinsic 
elongational  viscosity,  v,  defined  by 

v  =  (-q^)c  -►  0  ^  '2'9^ 

and,  for  elongated  ellipsoids,  is  depicted  in  Figure  (1.1)  as  a  func¬ 
tion  of  AR  with  a  as  a  parameter.  The  limiting  cases  defined  by 
a  =  0  and  a  =  °°  are  physically  attributable  to  particle  orientation 
in  an  extensional  flow  field.  In  the  case  of  a  =  0,  the  particles 
retain  a  completely  random  orientation  throughout  the  flow  field  as  a 


nsp  w  I/^dt' 
e  no^ 


where 


n  -  n. 


sp 
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FIGURE  (1.1)  ELONGATIONAL  VISCOSITY 

BEHAVIOUR  OF  ELONGATED 
ELLIPSOIDS 
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consequence  of  the  forces  of  Brownian  motion  being  large  in  comparison 
to  the  alignment  forces  of  the  flow  field.  In  the  case  of  a  =  ”, 
where  the  effect  of  brownian  motion  is  negligible,  the  particles  are 
oriented  parallel  to  the  streamlines.  This  is  a  condition  of  maximum 
energy  dissipation  (or  maximum  elongational  viscosity)  which  is 
approached  asymptotically  as  the  deformation  rate  tends  to  large 
values  for  a  suspension  of  very  small  particles. 

The  effect  of  brownian  motion  on  a  suspension  of  relatively 
large  particles,  such  as  fibre  suspensions,  will  be  negligible  except 
in  the  limit  of  zero  flow.  Hence,  the  intrinsic  elongational 
viscosity  of  a  suspension  of  large  particles  will  be  defined  by  the 
limiting  case  of  a  =  »  for  extensional  flows  of  interest.  This  sug¬ 
gests  that  the  elongational  viscosity  of  long  particles  will  be  very 
high,  but  essentially  free  of  elastic  behaviour. 

Determination  of  the  elongational  or  tensile  stress,  t^,  in 
purely  extensional  flows  of  dilute  suspensions  of  elongated,  rigid 
particles  is  given  by 

t-j-j  s  n  e -j i  (1 .2.10) 

and  independent  of  whether  or  not  the  suspension  exhibits  elastic 
behaviour.  It  is  significant  that  the  tensile  stresses  will  be  very 
large  for  long  particles.  Since  the  analysis  leading  to  Equation 
(1.2.10)  is  defined  only  for  purely  elongational  flow,  it  is 
necessarily  restricted  to  relating  the  total  elongational  stress  to 
the  elongational  deformation  rate  and  applicable  only  to  purely 
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extensional  flows  or  those  which  may  be  reasonably  approximated  as 
purely  extensional.  Solution  to  the  general  flow  problem  requires  a 
description  of  the  continuum  behaviour  through  the  use  of  a  constitu¬ 
tive  equation  which  relates  the  material  functions,  x..,  to  the 

'  J 

kinematic  variables,  d... 

'  J 

It  is  apparent  that  the  development  of  Takserman-Krozer  and 
Ziabicki  may  be  useful  in  defining  a  "better"  constitutive  relation. 
However,  since  this  study  is  concerned  with  the  behaviour  of  fibre 
suspensions  in  extensional  flow,  where  it  is  possible  to  relate 
elongational  stress  to  elongational  deformation  rate,  the  definition 
of  a  "better"  constitutive  relation  is  left  to  other  investigators. 
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CHAPTER  II 
THEORY 


2.1  GENERAL 

Unlike  measurements  of  shear  stresses,  the  measurement  of 
normal  stresses  involves  the  total  stress,  t.,-,  which  includes  the 

I  J 

isotropic  stress,  p6 . • ,  as  well  as  the  deviatoric  stress,  S • •  (27). 

■  vJ  *  J 

Since  the  meaning  of  p6^j  is  not  fundamentally  clear  for  visco¬ 
elastics  (24,25,27),  it  is  often  convenient  to  remove  the  influence 
of  p  from  normal  stress  measurements.  This  is  commonly  achieved  by 
presenting  results  in  terms  of  stress  differences,  for  example  (27): 

t-|i  ~  x 22  =  "  P  ^ii  “  (“P  S^)  ~  si -j  ■*  $22  (2.1.1) 

When  knowledge  of  the  total  stress  is  required,  it  is  necessary  to 
make  some  assumptions  about  p  (24,26,27). 

Although  measurements  have  been  made  of  normal  stresses 
developed  in  very  slow  elongational  flows  (3,24,25)  and,  more  recently, 
in  rapid  elongational  flows  (25),  most  measurements  have  involved 
normal  stresses  developed  in  simple  laminar  shear  flows  (hereinafter 
referred  to  as  SLSF).  Normal  stresses  developed  in  SLSF  of  poly¬ 
meric  solutions  can  be  measured  by  a  variety  of  techniques  (27). 

The  cone  and  plate  instrument  (26,27)  is  very  popular  due  to  the 
superior  flow  kinematics  and  to  the  absence  of  crude  assumptions 
necessary  to  complete  an  analysis  (27).  Another  popular  system  is 
the  capillary  jet  (26,27),  which  is  capable  of  producing  first  normal 
stress  data  at  relatively  high  shear  rates.  These  techniques,  however, 


. 
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are  practically  restricted  to  normal  stress  measurements  of  fairly 
elastic  polymeric  solutions  of  moderate  concentrations  (26,27). 

In  the  case  of  highly  concentrated  solutions  or  polymer  melts, 
for  example,  the  instabilities  which  develop  in  the  cone  and  plate 
geometry  limit  normal  stress  measurements  to  uninterestingly  low 
deformation  rates  (26).  Similarly,  the  capillary  jet  technique  is 
unwieldy  with  such  materials  and  subject  to  flows  of  low  Reynolds 
number  where  end  effect  corrections  may  be  appreciable  (2,14,16,17, 
26,28).  At  the  other  extreme  of  very  dilute  elastic  solutions,  or 
even  moderate  concentrations  of  slightly  elastic  solutions,  the 
low  levels  of  normal  stresses  developed  in  SLSF  become  difficult  to 
measure  with  any  degree  of  resolution  (26). 

Metzner,  Uebler,  and  Chan  Man  Fong  (26)  made  use  of  the  well- 
known  Maxwell  model  (27)  in  an  attempt  to  predict  normal  stresses  in 
the  entrance  region  of  long  circular  ducts.  Although  unsuccessful 
quantitatively,  they  were  able  to  suggest  that  sensitive  measurements 
of  normal  forces  could  be  realized  by  the  measurement  of  the  thrust, 
or,  equivalently,  the  expansion  of  orifice  jets.  This  suggestion 
was  made  on  the  basis  of  some  important  observations,  namely: 

i)  In  the  flow  of  elastic  fluids  through  a  contraction,  the 
kinematics,  immediately  upstream  of  the  contraction,  are 
essentially  those  of  a  purely  elongational  flow  (41). 
ii)  The  magnitude  of  the  normal  stresses  developed  in  these 
converging  flows  were  shown  to  be  large. 


'  *■'  Ki  1  ■ "  ’ 
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iii)  The  rate  of  elongational  deformation,  e,  in  these  converg¬ 
ing  flows  appeared  to  be  guided  by  the  constancy  of  the 
dimensionless  elongation  rate,  e^e  (37).  Here,  e.^ 
represents  the  relaxation  time  of  the  fluid  (35,36,37). 

The  orifice  jet  thrust  device  was  employed  by  Metzner  and  Metzner  (25) 
who  demonstrated  the  sensitivity  of  this  device  to  low  levels  of 
elasticity  as  well  as  the  high  tensile  stresses  developed  in  rapid 
elongational  flows  of  viscoelastic  media. 

2.2  ORIFICE  JET  THRUST  MEASUREMENTS 
2.2.1  Kinematics  of  Orifice  Flows 

Metzner  and  Metzner  (25)  have  defined  the  velocity  field  for  a 
viscoelastic  fluid  upstream  of  a  sharp-edged  circular  orifice  as 
depicted  in  Figure  (2.1).  The  only  upstream  fluid  which  passes 
through  the  orifice  is  confined  to  the  conical  region  bounded  by  the 
dashed  lines,  with  the  fluid  outside  this  conical  region  recirculating 
within  the  large  toroidal  vortices  indicated.  The  kinematics  of  this 
velocity  field  may  be  described  in  a  spherical  coordinate  system  hav¬ 
ing  its  origin  at  the  apex  of  the  cone  formed  by  converging  stream¬ 
lines.  In  this  coordinate  system,  the  physical  components  of  velocity 
may  be  approximated  by  the  equations  (25) 

V  =  - jr-S - 

r  2 it  r  (1  -  cos  $) 

v  =  v  =  0  (2.2.1 ) 

0 

where  Q  denotes  the  volumetric  flowrate,  4>  the  half-angle  of  the  cone 
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FIGURE  (2.1) 

VELOCITY  FIELD  FOR  VISCOELASTIC 
FLOW  THROUGH  AN  ORIFICE 
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and  r  the  distance  from  the  apex  of  the  cone. 

Equations  (2.2.1)  may  be  used  to  evaluate  the  physical  com¬ 


ponents  of  the  deformation  rate  tensor 

0 


d.  . 
ij 


9r 


0  r 
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(2.2.2) 


or 


d.  .  = 
13 
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irr  (l-COS$) 
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27ir  (1  -cos<J>) 
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2nr  ( 1 -cos$ ) 


(2.2.3) 


The  deformation  rates,  at  the  orifice,  can  now  be  defined  in  terms  of 
the  orifice  radius,  Rq: 


n  •  3 
-Q  sin  $ 

ttRq  ( 1  -COS4> ) 


3 

Q  sin  $ 

2nRg(l-cos<t) 


(2.2.4) 
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The  only  variable  in  Equations  (2.2.4)  which  may  be  difficult  to 
determine  experimental ly  is  $.  Since  the  dependence  of  $  on  Q  and 
Rq,  among  other  things,  is  unknown,  it  is  immediately  apparent  that 
$  cannot  be  measured  in  a  separate  (independent)  experiment  unless 
dimensional  similarity  is  preserved.  In  this  light,  the  observa¬ 
tions  of  previous  investigators  are  worth  consideration. 

Under  the  experimental  conditions  studied  by  Uebler  (41), 
the  full  cone  angle,  2$,  was  found  to  be  about  30°  and  independent 
of  flowrate.  More  recent  evidence  (24,25)  however,  suggests  that 
$  decreases  with: 

i)  increasing  flowrate, 

ii)  decreasing  orifice  size,  and 

iii)  increasing  elasticity  (an  indirect  observation  based  on 
the  constancy  of  the  dimensionless  elongational  rate, 

0fle)  (24,25,26). 

Uebler's  work  also  indicated  that  the  converging  flows  of  visco¬ 
elastic  media  were  essentially  free  of  shearing  deformations  (25,26,41) 
as  implied  by  Equations  (2.2.1).  Moreover,  Metzner  and  Metzner 
qualitatively  indicate  that  the  absence  of  shearing  deformations 
becomes  more  apparent  with  decreasing  $.  These  observations  show 
that  experimental  knowledge  of  4>  is  essential  to  accurately  predict 
deformation  rates.  Additionally,  in  the  orifice  jet  thrust  investiga¬ 
tion  of  fluids  of  a  low  level  of  elasticity,  the  constancy  of  e^e 
implies  that  $  will  be  large.  Consequently,  measurements  of  the 
velocity  field  may  be  required  to  accurately  predict  deformation 
rates  and  to  determine  whether  or  not  shearing  deformations  are 


. 
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properly  negligible. 

Available  observations  on  Newtonian  flow  through  a  contraction 
(41)  indicate  that  the  fluid  which  passes  through  the  contraction  is 
bounded  by  the  geometry  of  the  reservoir.  The  toroidal  vortices 
common  to  flows  of  viscoelastic  media  in  similar  geometries  are 
absent  in  the  Newtonian  case  (26).  When  the  contraction  is  a  sharp- 
edged  circular  orifice  in  a  flat  plate,  one  would  expect  the  upstream 
velocity  field  to  be  as  depicted  in  Figure  (2.2).  This  velocity  field 
may  be  approximated  (in  spherical  coordinates)  by  modifying  Equations 
(2.2.1)  to 


v  U)  =  - 5-^ -  •  f(<f>) 

2nr  (l-cos$) 


(2.2.5) 


where  f(<j>)  will  be  a  maximum  at  <p  =  0  and  identically  zero  at  <j>  =  90°. 
This  simply  indicates  that  the  flow  is  not  purely  elongational  and 
the  deformation  rate  tensor  will  have  non-zero  non-diagonal  elements. 
Since  <i>  =  90°, the  elongational  deformation  rates  will  be  larger  than 
those  for  viscoelastics.  Exact  numerical  solutions  to  velocity  fields 
upstream  of  a  contraction  are  available,  but  restricted  to  flows  of 
low  Reynolds  number  (9,42).  Although  Equations  (2.2.5)  are  a  rather 
poor  representation  of  the  velocity  field  depicted  in  Figure  (2.2), 
they  will  serve,  nevertheless,  to  define  a  maximum  conceivable 
deformation  rate  for  orifice  flows. 
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FIGURE  (2.2) 

VELOCITY  FIELD  FOR  NEWTONIAN 
FLOW  THROUGH  AN  ORIFICE 
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In  view  of  the  presence  of  elasticity  in  dilute  fibre  suspen¬ 
sions  (10,29),  one  would  expect  the  velocity  field  description  upstream 
of  an  orifice  to  be  bounded  by  the  extremes  describing  very  visco¬ 
elastic  and  Newtonian  orifice  flows.  Again,  the  implied  constancy 
of  e^-j e  suggests  that,  if  elasticity  is  present,  the  expected  deforma¬ 
tion  rates  for  fibre  suspensions  will  be  lower  than  those  of  Newtonians. 
In  any  event,  there  is  no  evidence  to  suggest  that  the  deformation 
rates  will  be  higher  than  those  of  Newtonians  in  orifice  flows.  Since 
no  data  are  available  to  describe  the  velocity  field  in  orifice 
flows  of  fibre  suspensions,  reliable  elongational  deformation  rates 
cannot  be  predicted.  The  need  for  velocity  measurements  is  clearly 
indicated. 


2.2.2  Analysis  of  Orifice  Jet  Thrust 

In  the  interpretation  of  the  stress  levels  of  the  flow  depicted 
in  Figure  (2.1),  it  is  convenient  to  employ  a  cylindrical  coordinate 
system,  r,  e,  z,  with  the  z-axis  directed  along  the  centerline  of  the 
velocity  field.  A  momentum  balance  between  a  section  far  upstream  of 

t 

the  orifice  and  the  fully-developed  jet  yields  the  following  equation 
for  thrust  (25) : 


T 


(p  V  2  -  T„)  d  A 

A 


(2.2.6) 


in  which  T  denotes  the  (measured)  total  thrust,  p  the  fluid  density, 
v^  the  axial  component  of  velocity  and  the  total  axial  stress  in 
the  fluid  at  the  orifice.  If  the  velocity  profile  is  flat  (cylindri¬ 
cal  coordinates),  then  t-j-|  is  independent  of  radial  position,  and 
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Equation  (2.2,6)  may  be  simplified  by  replacing  vz  by  VQ  and  removing 
the  integral,  giving 

T  -  Vo2  -TU>*Ro2  <2-2-7> 

where  VQ  is  the  average  velocity  at  the  orifice.  Since  the  radial 
stress,  t£2>  is  zero  in  the  wel 1 -developed  jet  (25),  Equation  (2.2.7) 
may  be  expressed  in  terms  of  stress  differences  (25)  as 

T11  -  S11  -  S22  “  P  Vo2  -  T/"Ro2 

=  PtQ-y)2--^  (2.2.8) 

7T  R  L  TtR/ 

0  0 

showing  the  measurement  of  thrust,  T,  as  a  function  of  flowrate,  Q, 
relates  the  normal  stress  difference,  -  S^,  to  the  elongational 
deformation  rate  defined  by  Equation  (2,2.4). 

If  n  is  required,  then  the  total  stress,  t^,  is  utilized 
rather  than  the  stress  difference.  Hence,  noting  that 

In  =  n  l  (2.2.9) 


and,  approximately,  that 

V  sin3  <t> 

•  _  o _ 

e  ”  R  (1  -  COS  4») 


(2.2.10) 


Equation  (2.2.7)  may  be  written  as  a  second-order  polynomial  in  VQ, 
namely 


sin3  $  v 

R  (1  -  cos  $)  vo 
ov 


n 


(2.2.11) 
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The  elongational  viscosity  of  viscoelastics  can  be  determined  from 

experimental  measurements  of  T,  V  ,  and  <t>.  The  form  of  Equation 

(2,2. 11)  suggests  a  simple  least-squares  approach  to  data  smoothing 

.  3 

provided  that  the  term  7TS'ln — ^  remains  essentially  constant  over 

U-cos 

the  range  of  an  experiment.  The  least-squares  approach  will  be  used, 
in  fact,  to  interpret  data. 

The  orifice  jet  thrust  of  an  ideal  flow  can  be  described  by  a 
momentum  balance  as  determined  by  Equation  (2.2.6).  In  this  case, 
the  axial  stress,  r-| -j ,  will  be  zero.  Noting  the  presence  of  a  vena 
contracta  a  short  distance  downstream  of  the  orifice,  the  thrust  will 
be  related  to  the  velocity  at  the  vena  contracta,  V^,  rather  than  on 
the  average  velocity  at  the  orifice,  V  .  Hence,  assuming  incompres¬ 
sible  flow,  the  thrust  is  given  by 


-  ^  Vo2  ■  “P  V  <2-2-12> 

where  is  the  cross-sectional  area  of  the  jet  at  the  vena  contracta. 
The  value  of  k  has  been  determined  from  ideal  flow  theory  to  be  1.64 
(5  ,22 ,25) * 

Since  orifice  flows  of  viscous  Newtonians  are  not  ideal, 
unusual  thrust  reduction  may  result  as  a  consequence  of  profile  relaxa¬ 
tion  and  surface  tension  effects  which  are  common  to  capillary  flows 
at  low  Reynolds  numbers  (2,14,16,17,28).  The  thrust  reduction  may 
also  be  affected  by  the  large  tensile  stresses  developed  immediately 
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upstream  of  the  orifice.  This  suggestion,  however,  hints  of  clair¬ 
voyance  since  Newtonian  liquids  have  no  memory  for  previous  deforma¬ 
tions  (24,25).  Hence,  the  thrust  per  unit  orifice  area  for  a  Newtonian 
orifice  jet  can  be  described  in  terms  of  the  orifice  velocity,  V  ,  as 

~TZ  kl  P  Vo2  '  k2  Vo  (2.2.13) 

1,1 Ro 

where  k-j  and  may  both  be  functions  of  velocity,  viscosity,  surface 
tension,  orifice  size  and  k^  may  also  be  a  function  of  deformation 
rate.  As  a  first  approximation,  it  is  assumed  that  k-j  and  k^  are 
constants,  independent  of  velocity. 

Noting  that  Equations  (2.2.11)  and  (2.2.13)  are  of  the  same 
form  (i.e.  both  in  terms  of  VQ  and  V  )  and  that  the  kinematics  of 
orifice  flows  of  fibre  suspensions  are  expected  to  lie  somewhere 
between  the  extremes  of  very  viscoelastic  and  purely  Newtonian,  the 
thrust  equation  for  suspension  flow  may  be  simply  written  without 
resorting  to  a  differential  momentum  balance.  Since  the  coefficient, 
k-j,  in  the  VQ  term  of  Equation  (2.2.13)  arises  from  the  non-flatness 
of  the  Newtonian  velocity  profile  at  the  orifice,  one  would  expect 
this  coefficient  to  be  smaller  for  fibre  suspensions  due  to  an 
expected  decrease  in  The  term  involving  VQ  in  Equation  (2.2.11) 
is  an  axial  stress  term  and  hence  necessary  in  the  fibrous  case. 
Similarly,  the  term  involving  VQ  in  Equation  (2.2.13)  is  a  consequence 
of  "unusual  behaviour"  (which  also  may  be  regarded  as  an  axial  stress) 
and  is,  therefore,  important  in  the  fibrous  case.  Thus,  the  orifice 
jet  thrust  of  a  fibrous  flow  will  be 
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(2.2.14) 


where  k-j  and  are  functions  similar  to  k-j  and  k^  respectively,  and 
k^  will  reflect  the  increase  in  the  axial  stress  (or  elongational 
viscosity)  due  to  the  presence  of  fibres. 

Although  Equations  (2.2.13)  and  (2.2.14)  both  appear  amenable 
to  a  polynomial  least-squares  regression  of  the  thrust  data,  a  "good" 
statistical  correlation  does  not  necessarily  imply  that  k-j  and  k^  (or 
k-j  and  k£)  represent  the  absolute  magnitude  of  the  actual  physical 
phenomena.  A  relative  trend  may  be  indicated  by  the  coefficients 
if  a  "physical  reason"  supports  the  constancy  and  numerical  magnitude 
of  k-j  (or  k-j)  over  the  range  of  an  experiment.  Hence,  in  order  to 
properly  compare  the  elongational  or  "elastic"  behaviour  of  a  fibre 
suspension  to  a  Newtonian,  it  is  necessary  to  compare  experimental 
results  for  a  Newtonian  and  a  suspension  which  have  approximately 
the  same  shear  viscosity  and,  possibly,  equal  surface  tensions. 
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CHAPTER  III 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 


3.1  GENERAL 

The  orifice  jet  thrust  technique,  noted  for  its  sensitivity  to 
low  levels  of  elasticity  (25,26),  was  chosen  for  this  study.  Sugar 
solutions  were  chosen  as  a  suspending  medium  for  the  fibres.  Since 
fibres  disperse  easily  in  a  medium  of  low  viscosity  and  particle 
aggregation  and  settling  are  minimal  in  a  medium  of  high  viscosity, 
suspension  preparation  was  simplified  by  first  dispersing  the  fibres 
in  water,  then  adding  sugar  to  increase  the  viscosity  of  the  suspend¬ 
ing  medi urn. 

The  experimental  program  consisted  of  measuring  thrust,  pres¬ 
sure  and  flowrate  for  orifice  flows  of  two  sugar  solutions  and  nine 
different  suspensions  of  fibres  in  sugar  solutions.  Originally, 
there  was  an  intention  to  make  velocity  measurements  upstream  of  the 
orifice  by  a  photographic  technique  similar  to  those  used  by  a  number 
of  others  (32,35,41).  However,  velocity  measurements  were  abandoned 
due  to  the  translucence  of  the  suspensions  used.  The  thrust  device 
and  ancillary  equipment  are  represented  schematically  in  Figure  (3.1) 
and,  in  some  detail,  in  Figures  (3.2),  (3.3)  and  (3.4). 

3.2  EXPERIMENTAL  APPARATUS 
3.2.1  Spring  System 

Basically,  the  system  consisted  of  the  orifice  device  supported 
by  two  32  inch  flat,  parallel  springs  clamped  firmly  at  each  end. 

The  spring  leaves  were  1/16  inch  thick  x  2  inches  wide  304  stainless 
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SCHEMATIC  OF  EXPERIMENTAL  APPARATUS 
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FIGURE  (3.2) 

TOP  VIEW  OF  THRUST  DEVICE 
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SIDE  VIEW  OF  UNDERCARRIAGE  OF  SPRING  TABLE 
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steel,  and,  when  clamped  to  the  supports,  could  leave  a  maximum  of 
31  inches  of  undamped  length.  The  two  5-3/4  inches  wide  by  6-3/4 
inches  high  rectangular  spring  supports,  shown  in  Figure  (3.2), 
were  milled  from  a  1  inch  cold  roll  mild  steel  plate.  The  spring 
supports  were  attached  to  the  1  inch  x  1  inch  x  34  inch  mild  steel 
support  bar  by  means  of  the  spring  support  clamps,  which  were  welded 
to  the  spring  supports.  Each  clamp  was  fabricated  from  a  2  inch 
cube  of  mild  steel  and  fixed  rigidly  to  the  support  bar  by  two  No. 

10  -  24  machine  set  screws  as  shown  in  Figure  (3.2).  The  support 
bar  was  bolted  rigidly  to  a  20  inch  x  40  inch  x  3/4  inch  thick 
rectangular  aluminum  base.  Each  side  of  the  spring  supports  was 
drilled  and  tapped  to  accept  two  No.  10  -  24  machine  bolts  on  2-1/2 
inch  centres  with  the  top  hole  located  1/2  inch  below  the  top  of  the 
support.  This  facilitated  clamping  the  spring  leaves  to  the  sides 
of  the  supports  by  means  of  four  1  inch  x  3-1/2  inch  x  1/4  inch  flat 
mild  steel  plates  which  were  bolted  to  the  supports. 

A  small  table  of  mild  steel,  1/4  inch  x  3  inches  x  7  inches 
in  size,  was  mounted  above  the  springs,  halfway  along  their  unclamped 
length,  with  the  long  axis  of  the  table  perpendicular  to  the  spring 
face.  As  shown  in  Figure  (3.3),  the  spring  table  was  attached  to 
the  springs  by  means  of  spring  table  supports.  Each  of  the  two  spring 
table  supports  consisted  of  two  lengths  of  3/4  inch  half-rounds  whose 
rounded  surfaces  clamped  a  spring  leaf  between  them.  The  longer, 

3-1/4  inch,  half-rounds  were  rigidly  bolted  to  the  spring  table. 

The  viscous  damper  depicted  in  Figure  (3.3)  consisted  of  a 
round  paddle,  which  was  rigidly  fixed  to  the  spring  table  supports  by 
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means  of  a  3/4  inch  x  3/16  inch  thick  mild  steel  strap,  and  immersed 
in  a  trough  of  S.A.E.  140  gear  oil.  The  trough  was  a  fabricated 
4  inch  x  1  inch  I.D.  cylinder  with  enclosed  ends  and  a  cutout  top  to 
accommodate  the  paddle.  The  trough  was  rigidly  attached  to  the 
aluminum  base  by  means  of  1-1/4  inch  x  1/4  inch  mild  steel  strap. 

The  entire  apparatus  was  enclosed  by  a  19  inch  x  39  inch  x  18 
inch  x  1/4  inch  thick  plexiglass  box  with  a  removable  top.  Small 
openings  were  provided  for  the  transducer  leads  and  the  feed  tubing. 

A  large  hole  (2-1/2  inch  <j>)  was  cut  opposite  the  orifice  plate  to 
accommodate  a  short  length  of  2  inch  I.D.  Tygon  tubing  which  served 
to  collect  the  fluid  leaving  the  orifice. 

3.2.2  Orifice  Plate  and  Reservoir 

i 

The  orifice  reservoir,  which  was  bolted  to  the  spring  table, 
was  fabricated  from  3/8  inch  plexiglass.  Sheets  of  the  material  were 
glued  and  screwed  together  to  form  a  4  inch  x  4  inch  x  8  inch  open 
box.  The  orifice  plate  was  fabricated  from  1/4  inch  brass  plate  and 
bolted  to  one  end  of  the  open  box.  A  sharp-edged  orifice,  bevelled 
30°  normal  to  the  plate,  was  located  in  the  centre  of  the  plate.  A 
feed  plate,  also  fabricated  from  1/4  inch  brass  plate,  served  as  the 
other  end  of  the  open  plexiglass  box.  The  feedplate  originally  held 
a  thin,  brass,  perforated  flow-distributor  plate  2  inches  into  the 
box  by  means  of  four  long  bolts.  Owing  to  fibres  being  trapped  in  the 
perforations ,  this  was  later  abandoned  in  favour  of  a  2  inch  x  2  inch 
x  1/4  inch  thick  unperforated  teflon  plate.  The  centre  of  the  feed 
plate  was  drilled  and  tapped  to  accept  the  0-ring  side  of  a  5/8  inch 
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O-ring  seal  x  3/8  inch  swagelok  fitting.  The  swagelok  fitting  held 
a  fabricated  3/8  inch  copper  tubing  elbow  to  which  was  clamped  the 
3/8  inch  heavy-wall  Tygon  feed  tubing. 

3.2.3  Transducers  and  Mountings 

One  side  of  the  orifice  reservoir  was  drilled  and  tapped  to 
accept  the  0-ring  side  of  a  3/8  inch  0-ring  seal  by  1/4  inch  swagelok 
fitting.  Instrument  tubing  connected  the  swagelok  fitting  with  the 
Validyne  DP75  differential  pressure  transducer,  manufactured  by 
Validyne  Engineering  Corporation  of  Northridge,  California.  The 
negative  side  of  the  transducer  was  open  to  the  atmosphere  while  the 
positive  side  was  connected  to  the  tee  and  valve  arrangement  depicted 
schematically  in  Figure  (3.4).  The  rigid  signal  line  support,  shown 
in  Figure  (3.2),  rigidly  fixed  the  instrument  tubing  to  the  spring 
support  and  served  to  isolate  accidental  movements  of  the  externally 
located  transducer  from  the  calibrated  spring  system.  The  calibration 
vessel,  which  was  connected  to  the  transducer  with  instrument  tubing, 
served  a  two-fold  purpose;  namely,  to  calibrate  the  transducer  (as 
discussed  later)  and  to  bleed  the  transducer  and  signal  lines  prior 
to  an  experimental  run. 

Spring  displacement  was  measured  by  a  Boulton  Paul  model  F51M 
zero-resistance  displacement  transducer,  manufactured  by  Boulton  Paul 
Aircraft  Ltd.  of  Wolverhampton,  England.  The  armature  of  the  trans¬ 
ducer  was  attached  to  the  rear  spring  table  support,  normal  to  the 
plane  of  the  spring  leaves.  The  transducer  coil  was  held  by  a  1-1/4 
inch  x  1-1/4  inch  x  1-1/2  inch  brass  prism  which  was  bored  out  and 
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slit  along  the  bore  to  serve  as  a  clamp  for  the  cylindrical  trans¬ 
ducer  housing.  As  shown  in  Figure  (3.3),  the  transducer  clamp  was 
rigidly  fixed  to  the  aluminum  base  via  1-1/4  inch  x  1/4  inch  mild 
steel  strap. 

Both  transducers  were  excited  and  demodulated  by  Hewlett- 
Packard  Sanborn  8805A  Carrier  Preamps.  The  circuit  of  the  pressure 
transducer  required  an  RC  circuit  to  correct  the  phase  shift  between 
input  and  output  signals.  The  displacement  transducer  required  no 
signal  conditioning.  Outputs  of  the  preamps  were  monitored  by  a 
two-channel  Dynagraph  type  RS  strip-chart  recorder  manufactured  by 
Beckman  Instruments. 

3.2.4  Pumps 

The  original  pumping  system  used  in  this  study  was  a  pair  of 
Zenith  type  5,  model  B-4391  gear  metering  pumps  with  teflon  seals, 
connected  in  parallel.  Power  was  supplied  to  the  pumps  by  a  1-hp. 
220V,  3600  rpm.  synchronous  motor  manufactured  by  Gill  Electric  Motors 
(GEL)  of  Brighton,  England.  The  pump  speed  was  regulated  by  a  60- 
speed  epicentric  transmission  manufactured  by  Sangamo  Controls  Ltd. 
of  Bognor  Regis,  England.  Although  this  system  was  able  to  deliver 
20  cc/sec.,  which  was  sufficient  for  this  study,  a  problem  with  fibre 
plugging  and  the  eventual  availability  of  an  Instron  tester  made  way 
for  a  more  convenient  pumping  system. 

The  latter  pumping  in  this  study  was  supplied  by  a  single- 
acting,  6  inch  I.D.  x  20  inch  stroke  hydraulic  cylinder  powered  by  an 
Instron  model  TTBM  tester.  This  system  provided  the  ideal  conditions 
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of  steady,  positive  displacement  and  pre-cal ibrated  flow. 

3.2.5  Experimental  Fluids 

The  fluids  used  in  this  study  are  classified  in  Table  (3.1) 
along  with  the  corresponding  test  numbers.  The  suspending  medium 
(for  the  fibres)  was  a  solution  of  Alberta  beet  sugar  in  water. 

A  60%  (by  weight)  sugar  solution  served  as  the  suspending 
medium  for  tests  2  to  5,  where  the  short  fibres  were  used.  However, 
an  attempt  to  obtain  thrust  measurements  for  the  longest  fibres 
suspended  in  the  60%  solution  was  aborted  due  to  fibre  aggregation 
in  the  orifice  reservoir.  It  was  found,  by  trial  and  error,  that 
a  66  2/3%  sugar  solution  would  suspend  even  the  longest  fibres  used 
in  this  study  without  noticeable  fibre  aggregation.  Hence,  a 
66  2/3%  sugar  solution  suspended  the  fibres  in  the  latter  tests 
(tests  8  to  16) , 

The  viscosity  of  the  66  2/3%  sugar  solution  was  sampled  and 
found  to  be  1.7  poise  at  22  ±  1°C.  This  agrees  quite  favourably 
with  the  handbook  value  (43)  of  1.8  poise  at  22°C.  The  shear 
viscosities  of  the  suspensions,  however,  were  not  measured  since 
the  effect  of  shear  viscosity  was  originally  considered  to  be  of 
minor  importance.  Moreover,  at  the  concentration  levels  of  this 
study,  the  change  in  shear  viscosity,  due  to  the  presence  of  fibres, 
was  small  (less  than  50%)  (29)  compared  to  the  150%  change  between 
the  60%  and  66  2/3%  sugar  solutions  which  resulted  in  an  almost 
negligible  difference  in  thrust  reduction  (see  Figure  (4.2)). 
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TABLE  (3.1) 
EXPERIMENTAL  FLUIDS 


Suspending  Fibre 

Fluid*  Axis  Ratio 

Fibre 

Concentrati  on 
(Wt.%)  (Vol.%) 

Measured 

Density 

(qm./cc. 

A 

- 

0.00 

0.00 

1 .28 

A 

85 

0.05 

0.045 

? 

A 

85 

0.10 

0.090' 

? 

• 

A 

85 

0.15 

0.135 

? 

A 

85 

0.20 

0.180 

? 

A 

- 

0.00 

0.00 

? 

• 

C 

- 

0.00 

0.00 

1.00 

B 

85 

0.30 

0.278 

1  .32 

B 

170 

0.10 

0.0926 

1  .32 

B 

170 

0.20 

0.185 

1  .32 

B 

255 

0.10 

0.0926 

1  .32 

B 

- 

0.00 

0.000 

1.32 

B 

340 

0.10 

0.0926 

1  .32 

B 

85 

0.20 

0.185 

1.32 

B 

- 

0.00 

0.000 

1.32 

B 

170 

0.10 

0.0926 

1.32 

*  Suspending 

Fluid  A: 

60%  Sugar  Solution 

Suspending 

Fluid  B: 

66-2/3%  Sugar  Solution 

Suspending 

Fluid  C: 

0.10%  Separan  AP30 

Sol ution 

ORIFICE  DIAMETER 

Tests 

1  to  4  = 

0.0846 

cm. 

Tests 

5  to  16  = 

0.1477 

cm. 

"  a 
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As  indicated  by  Table  (3d),  suspension  densities  were  not 
measured  in  tests  2  to  6.  At  this  point  in  the  study,  it  was  found 
that  evaporation,  during  suspension  preparation  over  two  or  three 
days  at  high  temperature,  caused  a  water  loss  of  up  to  40%.  The 
amount  of  water  lost,  which  was  unknown  in  tests  2  to  6,  represented 
a  change  in  suspension  density,  which,  in  turn,  classified  these 
tests  as  unreliable,,  Since  the  thrust  reduction,  due  to  the  pre¬ 
sence  of  fibres,  was  relatively  small  (see  Chapter  IV),  it  is 
emphasized  that  a  small  (and  unknown)  change  in  suspension  density 
will  significantly  affect  the  absolute  thrust  reduction.  In  tests 
8  to  16,  the  suspension  densities  were  adjusted  (by  adding  water) 
to  the  required  specification  (see  comments  in  following  sections). 

The  fibres  used  in  this  study  were  purchased,  precut,  from 
Microfibres  Inc.  of  Pawtucket,  Rhode  Island.  All  fibres  were  1.5 
denier  viscose  rayon  cut  to  lengths  of  0.05  inch,  0.10  inch,  0.15 
inch,  and  0.20  inch.  Sample  measurements  of  dry  fibre  diameter, 
with  the  aid  of  a  measuring  microscope,  yielded  a  remarkably  con¬ 
stant  value  of  15  ±  0.25  microns.  Examination  of  the  fibres,  with 
a  microscope,  after  a  run,  revealed  no  evidence  of  frilibration  or 
swelling.  Vernier  caliper  measurements  of  a  sample  of  0.05  inch 
fibres  showed  the  stated  lengths  to  be  accurate  within  limits  of 
this  measurement.  The  density  of  the  fibres  was  determined  by 
immersing  fibres  in  various  measured  mixtures  of  carbon  tetrachloride 
and  alcohol  until  no  gravitational  gradient  was  observed.  The 
density  thus  determined  was  found  to  be  1.425  gm/cc. 
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3.3  EXPERIMENTAL  PROCEDURE 

3.3.1  Suspension  Preparation 

Sugar  solutions  were  prepared  by  weighing  the  ingredients, 

and  stirring  and  heating  the  solution  to  expel  the  suspended  air 

bubbles.  The  solution  was  then  cooled  while  mixing  continuously. 

When  the  solution  cooled  to  room  temperature,  it  was  weighed  to 

determine  the  evaporation  loss  and  sufficient  water  was  added  to 

* 

make  up  this  loss  . 

In  preparation  of  fibre  suspensions,  the  required,  weighed 

amounts  of  fibres  and  water  were  first  vigorously  mixed  to  disperse 

the  individual  particles.  The  sugar  was  added  after  dispersion  was 

assumed  to  be  complete.  The  heating,  cooling  and  checking  procedure 

* 

was  then  followed  as  outlined  above  . 

The  partially  degraded  polymer  solution  used  in  test  7  of 
this  study  was  obtained  as  a  sample  from  a  large  batch  which  had 
been  prepared  for  pipeline  tests  (32). 

3.3.2  Calibration 

Flow  calibration  was  a  relatively  simple  procedure  for  both 
pumping  systems,  since  the  flowrates  could  be  easily  related  to 
known  pump  speeds.  The  calibration  was  accomplished  by  weighing 
the  fluid  collected  over  a  period  of  about  2  minutes.  Flow  calibra¬ 
tion  data  are  presented  in  Appendix  A. 

Calibration  of  the  pressure  transducer  was  accomplished 


Evaporation  losses  were  adjusted  only  in  tests  8  to  16. 


26 

q  ns9d 
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through  the  use  of  the  scheme  represented  in  Figure  (3.4).  The  cali¬ 
bration  vessel  was  a  cylindrical  bomb  (closed  at  both  ends)  which  was 
half-filled  with  water.  An  opening  in  the  top  of  the  bomb  could  be 
connected  to  either  a  mercury  manometer  or  a  water  manometer,  depend¬ 
ing  upon  the  desired  calibration  range.  A  line  from  the  bottom  of 
the  vessel  to  valve,  V2,  provided  the  pressure  signal  to  the  transducer. 
Prior  to  a  calibration,  valve  VI,  was  closed,  V2  was  opened,  and  the 
transducer  bled  and  zeroed.  Calibration  was  performed  by  pressurizing 
the  bomb  to  various  levels  and  reading  a  manometer  and  corresponding 
signal  on  the  recorder.  Calibration  data  are  presented  in  tabular 
and  graphic  form  in  Appendix  A. 

The  thrust  calibration  was  accomplished  by  attaching  one  end  of 
a  thin  cotton  thread  to  the  rear  spring  table  support  and  the  other 
end  to  an  external  support  located  behind  and  above  the  thrust 
apparatus  in  a  plane  perpendicular  to  the  axis  of  the  springs.  Balance 
weights  were  attached,  by  a  short  length  of  thread,  to  this  thread, 
near  the  thrust  apparatus  and  an  arbitrary  weight  attached  by  a  length 
of  thread  to  the  original  thread,  farther  away  from  the  thrust  device. 
The  arbitrary  weight  was  immersed  in  a  container  of  water  to  damp 
oscillations.  The  original  thread  was  adjusted  such  that  the  thread 
coming  from  the  spring  table  support  was  level  (perpendicular  to  the 
plane  of  the  springs).  A  measurement  of  the  hypotenuse  formed  and  the 
horizontal  distance  between  the  known  and  arbitrary  weights  enabled 
solving  the  "force-triangle"  for  the  horizontal  force  on  the  springs. 

A  30  cm.  stainless  steel  machinist  rule,  with  1  mm.  graduations,  was 
used  to  measure  the  hypotenuse  and  horizontal  distances. 
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The  force  calibrations  were  used  directly  instead  of  relating 
them  to  spring  displacement.  Since  the  spring  calibration  was  sensi¬ 
tive  to  the  way  the  pressure  and  feed  lines  were  connected,  it  varied 
from  test  to  test.  Hence,  a  spring  calibration  was  performed  for  each 
test.  Calibration  data,  their  least-squares  equations,  and  a  brief 
discussion  on  the  errors  involved  are  presented  in  Appendix  A. 

3.3.3  Data  Collection 

The  equipment  was  assembled,  pressure  transducer  bled,  and  both 
transducers  zeroed.  An  experimental  test  was  performed  through  the 
flowrates  as  given  in  the  data  (Appendix  B).  The  pressure  and  thrust 
recorded  for  each  flowrate,  as  well  as  the  orifice  size,  are  also  pre¬ 
sented  in  Appendix  B.  Upon  completion  of  a  test,  the  thrust  calibration 
was  performed.  All  tests  were  performed  at  a  temperature  of  22°  1  1°C. 

Orifice  "plugging"  was  observed  when  the  orifice  jet  ceased 
and  the  pressure  signal  rose  monotonical ly  (to  the  point,  often,  where 
the  plexiglass  reservoir  burst).  In  cases  where  the  reservoir  did  not 
burst,  the  "plugging"  problem  was  relieved  be  carefully  reaming  the 
orifice  with  a  short  length  of  thin  wire.  Tests  where  "plugging" 
occurred  were  invariably  unsuccessful  as  a  result  of  either  a  broken 
reservoir,  or  a  calibration  change  in  the  spring  system.  Hence,  all 
tests  where  "plugging"  occurred  were  deemed  unreliable. 
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CHAPTER  IV 

DISCUSSION  OF  RESULTS 

4.1  VISCOELASTIC  BEHAVIOUR 

Thrust  measurements  of  the  0.10%  Separan  AP30  solution  are  pre¬ 
sented  here  to  demonstrate  the  sensitivity  of  the  orifice  jet  thrust 
device  and  to  compare  with  previous  results.  Figure  (4.1)  depicts  the 
results  of  test  7  as  well  as  the  results  of  Metzner  and  Metzner  (25) 
for  water  and  0.083%  Separan  AP30  solution.  The  results  of  test  13 
are  also  presented  to  show  the  relatively  low  level  of  "elasticity" 
developed  in  the  most  "elastic"  suspension  used  in  this  study.  The 
results  of  test  7  are  generally  in  good  agreement  with  Metzner  and 
Metzner  for  a  similar  concentration  of  polymer. 

4.2  NEWTONIAN  BEHAVIOUR 

Anticipation  of  unusual  thrust  reduction  (or,  effectively  non¬ 
ideal  flow)  in  orifice  flows  of  viscous  Newtonians  motivated  orifice 
jet  thrust  measurements  for  a  Newtonian  fluid  of  approximately  the 
same  shear  viscosity  as  the  fibre  suspensions.  Thus,  it  was  possible 
to  determine  any  change  in  elongational  viscosity  resulting  from  the 
fibres.  As  shown  in  Figure  (4.2),  the  thrust  per  unit  area  of  the 
Newtonian  suspending  medium  deviated  to  lower  values,  at  lower 
velocities,  than  those  of  water  (25).  It  is  known  that  surface  ten¬ 
sion  becomes  important  in  viscous  capillary  flows  at  low  velocities 
and  tends  to  significantly  reduce  thrust  (2,14,16,17,28).  Surface 
tension  may  also  tend  to  reduce  the  thrust  in  viscous  orifice  flows 
in  addition  to  the  thrust  reduction  caused  by  the  axial  stress  (see 
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VELOCITY  (cm./sec.) 


FIGURE  (4.1) 

THRUST- VELOCITY  CURVES  FOR  WATER, 
SEPARAN  AP30  SOLUTIONS  AND  THE 
MOST  "ELASTIC"  SUSPENSION 
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THRUST- VELOCITY  CURVES  FOR  WATER 
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Section  (2.2)).  The  data  presented  in  Figure  (4.2)  do  not,  however, 
distinguish  the  singular  importance  of  surface  tension  from  the  overall 
thrust  reduction  in  orifice  jets. 

Thrust  measurements  of  viscous  Newtonian  flow  through  an  orifice 
are  presented  here  to  serve  as  a  comparison  for  thrust  measurements  of 
fibre  suspensions.  Since  the  difference  in  thrust  between  the  Newtonian 
and  suspension  flows  was  small,  they  will  not  be  presented  in  the  form 
used  by  Metzner  and  Metzner.  Instead,  they  will  be  normalized  to  an 
"averaged"  Newtonian  thrust  and  presented  on  the  basis  of  a  normal¬ 
ized  thrust,  TN.  However,  before  such  a  presentation  can  be  made, 
some  averaging  criterion  must  be  established. 

The  form  of  Equation  (2.2.13)  immediately  suggests  the  use  of  a 
polynomial  least-squares  regression  to  fit  the  flowrate  -  thrust  data. 

It  was  found  that  a  Type  III  fit  (see  Appendix  C)  best  represented 
the  Newtonian  thrust  data.  Hence,  the  average  thrust  per  unit  orifice 
area  of  a  Newtonian  is  given  by: 


(-$  -  CM  C  -  [k,]  v0 

ttK7 


(4.2.1) 


NEWT 


where 

< 


represents  the  averaged  Newtonian  thrust  per  unit 


NEWT 


orifice  area  and  [k-j]  and  [k2]  are  least-squares  estimates  of  the 
coefficients  k^  and  k2  defined  by  Equation  (2.2.14).  The  normalized 
thrust,  for  a  given  data  point,  will  then  be  defined  by 
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(4.2.2) 


An  averaged  normalized  thrust,  for  a  given  test,  can  be  defined  as 


Ctc-j]  Vq  -  [4]  v£ 


(4.2.3) 


Where  (tn)  is  simply  a  comparison  of  the  smoothed  thrust  data,  for  a 
given  test,  to  the  best-fit  curve  of  the  Newtonian  data.  The  para¬ 
meters  [k-j]  and  [4]  were  determined  by  a  Type  III  fit  for  each  test 
performed  in  this  study.  The  Newtonian  coefficients  [4]  and  [4] 
were  determined  by  a  Type  III  fit  on  the  combination  of  tests  1  and 
6  for  a  pure  60%  sugar  solution  (see  Table  3.1)).  It  was  found  that 
the  thrust  per  unit  orifice  area  could  be  expressed  as 


=  1  .767  VQ2  -  163.9  VQ  (4.2.4) 

nR0 

2 

where/ — ^)is  given  in  dyne/cm;  and  VQ  in  cm. /sec.  Similarily,  the 

ttR 

results  ot  the  combination  of  tests  12  and  15  showed  that  the  thrust 
per  unit  orifice  area  for  the  66-2/3%  sugar  solution  could  be 
expressed  as 

L)  =  1.8667  VQ2  -  220.14  VQ  (4.2.5) 

*Ro 

Figure  (4.3)  represents  TN  and (TN) for  tests  12  and  15.  The  increased 
scatter  in  TN  for  lower  velocities  suggest  that  experimental 
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FIGURE  (4.3) 

NORMALIZED  THRUST  FOR  TESTS  12  AND  15 
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uncertainty  increases  at  lower  velocities.  This  uncertainty  was  anti¬ 
cipated  since  the  temperature  sensitivity  of  the  spring  system  made 
small  measurements  especially  difficult.  The  apparent  scatter  in  TN 
at  high  velocities,  though  regarded  here  as  experimental  uncertainty, 
can  be  qualified,  in  part,  by  the  existence  of  an  unusual  behaviour 
which  will  be  discussed  in  Section  (4.4).  Although  the  reproducibility 
of  Newtonian  data,  at  high  velocities,  appears  to  be  within  1%  on  the 
basis  of  ^TN^,  the  scatter  in  TN  is  a  less  biased  indicator  of  the 
reproducibility.  Nevertheless,  the  close  agreement  of  ^TN)  with 
unity  demonstrates  the  utility  of  the  least-squares  smoothed  data 
in  the  presentation  of  results. 

4.3  FIBRE  SUSPENSION  BEHAVIOUR 
4.3.1  Thrust  Measurements 

In  principle,  the  normal  stress  difference,  s i q ”^22 9  ™  ^bre 
suspension  flows  can  be  calculated  from  Equations  (2.2.8)  and  may  be 
representative  of  an  increase  in  elongational  viscosity  at  high  flow- 
rates.  Low  flowrate  data  is  essentially  useless  since  experimental 
uncertainty  masks  the  trends  of  the  absolute  thrust  measurements.  At 
high  flowrates,  the  normal  stress  calculation  is  reliable  only  in  a 
relative  sense  due  to  the  uncertainty  in  the  magnitude  of  the 
coefficient  k-j  in  Equation  (2.2.14).  In  addition,  the  absence  of 
deformation  rate  data  leaves  no  convenient  variable  with  which  to 
correlate  the  normal  stresses.  Hence,  the  thrust  data  alone  will 
qualitatively  describe  the  increase  in  normal  stress  (or,  equiva¬ 
lently,  the  elongational  viscosity). 


. 


. 
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As  noted  in  Section  (4.2),  the  relatively  small  differences  in 
thrust  between  Newtonians  and  fibre  suspensions  can  be  conveniently 
presented  in  terms  of  a  normalized  thrust,  TN.  Since  tests  2  to  6 
were  considered  unreliable  due  to  the  uncertainty  in  sugar  concentra¬ 
tion,  they  will  not  be  considered  here  except  for  test  3  which  will 
serve  as  a  comparison.  The  results  of  these  tests  have  been  normal¬ 
ized  and  appear,  in  normal  form,  in  Appendix  B.  Test  8  was  unsuc¬ 
cessful  due  to  plugging.  It  was  noted  in  Chapter  III  that  the  plug¬ 
ging  itself  was  not  a  serious  problem,  however,  a  calibration  change 
in  the  spring  system  resulted  from  the  unplugging  procedure.  Tests 
9  to  16  are  considered  reliable  and  will  be  discussed  further. 

Tests  9  to  16  and  their  least-squares  representations  were 
normalized  with  respect  to  Equation  (4.2.5).  Figures  (4.4)  to  (4.9) 
depict  TN  and  ^TN^  for  tests  9  and  16,  10,  11,  13,  14  and  3  respect¬ 
ively.  As  in  the  case  of  the  normalized  Newtonian  representation  of 
Figure  (4.3),  the  fibre  results  indicate  increased  scatter  in  TN  at 
lower  velocities.  The  scatter  is  attributable  to  experimental 
uncertainty  and  to  "partial -pi ugging"  which  was  apparent  only  at  the 
lower  velocities.  The  "partial -pi ugging"  became  especially  apparent 
for  the  longer  fibres  when  the  jet  would  not  issue  perpendicular  to 
the  orifice  plate.  This  effect  immediately  disqualifies  any  correla¬ 
tion  between  "elastic"  effects  and  behaviour  of  TN  or  ^TN^  at  the 
lower  velocities.  The  apparent  scatter  in  TN  at  higher  velocities 
can,  again,  be  qualified,  in  part,  by  an  unusual  behaviour  which 


will  be  discussed  later. 
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FIGURE  (4.4) 

NORMALIZED  THRUST  FOR  TESTS  9  AND  16 


TN  &  <TN>  (dimensionless) 
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FIGURE  (4.5) 

NORMALIZED  THRUST  FOR  TEST  10 
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<TN>  (dimensionless) 
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FIGURE  (4.6) 

NORMALIZED  THRUST  FOR  TEST  11 


<TN>  (dimensionless) 
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FIGURE  (4.7) 

NORMALIZED  THRUST  FOR  TEST  13 


<TN>  (dimensionless) 
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FIGURE  (4.8) 

NORMALIZED  THRUST  FOR  TEST  14 


TN  &  <TN>  (dimensionless) 
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FIGURE  (4.9) 

NORMALIZED  THRUST  FOR  TEST  3 
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Qualitatively,  the  data  show  a  decrease  in  thrust,  or,  equi¬ 
valently,  increased  normal  stresses,  with  increasing  fibre  concentra¬ 
tion  and  axis  ratio,  AR.  The  effect  of  AR  on  "elasticity"  is  depicted 
in  Figure  (4*10)  with  averaged  normalized  thrust  curves  from  previous 
figures*  Any  bias  inherent  in  the  averaging  criteria  is  of  little 
consequence  since  only  relative  arguments  are  used.  In  other  words, 
the  bias  (if  present)  will  affect  all  the  results  essentially  the 
same  way  and  since  only  relative  trends  are  expected,  they  will  be 
preserved* 

Except  for  test  16,  which  appears  to  exhibit  an  unusually  low 
thrust,  the  effect  of  increasing  axis  ratio  in  a  regular  fashion  is 
to  reduce  thrust  in  a  regular  fashion.  The  irregularity  in  test  16 
is  attributable  to  a  consistent  "unusually  low  behaviour".  This 
irregularity  is  discussed,  in  more  detail,  in  Section  (4.4), 

The  effect  of  fibre  concentration  on  {tn)  is  depicted  in 
Figure  (4*11)  by  a  comparison  of  test  3  with  test  14  and  a  comparison 
of  tests  9  and  16  with  test  10*  Although  the  concentration  effect  is 
not  as  pronounced  as  the  effect  of  axis  ratio,  it  is  qualitatively 
consistent  with  the  theory  of  specific  viscosity  (39,40).  In 
principle,  Equation  (2.2.14)  defines  the  added  axial  stress  (due  to 
fibres)  in  orifice  flows  of  suspensions.  At  present,  it  is  not  known 
whether  the  added  stress  is  viscous  or  elastic  in  nature.  This  added 
stress,  nevertheless,  represents  an  increase  in  the  elongational 
viscosity  which,  in  terms  of  shear  viscosity,  is  greater  than  the 
Newtonian  value  of  3y.  Estimates  of  elongational  viscosities  of  the 
fibre  suspensions  are  presented  in  the  following  section. 
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FIGURE  (4.10) 

EFFECT  OF  AR  ON  <TN> 
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FIGURE  (4.11) 

EFFECT  OF  C  ON  <TN> 
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A  comparison  of  the  contraction  coefficients  for  the  fibre  sus¬ 
pensions  with  the  Newtonian  suspending  medium  is  depicted  in  Figure 
(4.12).  The  comparison  merely  indicates  that  the  pressure  drop  for  a 
given  flowrate  of  suspension  is  approximately  the  same  as  the  Newtonian 
value  for  the  high  velocities  where  the  data  are  considered  reliable. 

At  lower  velocities,  where  "partial-plugging"  may  have  been  present, 
the  pressure  drop  would  be  too  high  (as  indicated  in  Figure  (4.12)), 
due  to  an  effective  reduction  in  orifice  area. 

4.3.2  Prediction  of  Elongational  Viscosity 

The  absence  of  deformation  rate  data  immediately  disqualifies 
any  attempt  to  accurately  predict  elongational  viscosities  of  the 
fibre  suspensions.  It  is  possible,  however,  to  crudely  show  the 
qualitative  effects  of  axis  ratio  and  fibre  concentration  on  the 
elongational  viscosity  and  to  compare  these  effects  to  those  pre¬ 
dicted  by  Takserman-Krozer  and  Ziabicki  (39,40).  This  can  be  done 
by  numerically  manipulating  the  least-squares  coefficients  of  the 
polynomial  representation  of  thrust. 

As  noted  in  the  previous  section,  the  thrust  per  unit  orifice 
area  of  a  viscous  Newtonian  can  be  adequately  represented  by 

J-2  -  [kp  V02  -  [k2]  V0  (4.3.6) 

"Ro 

The  equivalent  thrust  equation  for  a  suspension  of  fibres  is  given 

by 
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“S'  =  [ki]  vo2  •  Ck2]  v0  (4.3.7) 

"V  10  20 

Assuming  that  the  term  [k^] VQ  is  truly  indicative  of  the  tensile 
stresses  developed  in  the  extensional  flow  field,  a  "representative" 
axial  stress  is  given  by 

{t1 1 }  =  [k2]  VQ  (4.3.8) 

Since  n  =  3y  for  three-dimensional  Newtonians  (27),  a  "representative" 
deformation  rate,  { e } ,  can  be  defined  as 


(e)  = 


h)ll 

3y 


k2  Vo 

3u 


(4.3.9) 


The  term  "representative"  refers  to  numerically  relative  quantities 
as  opposed  to  actual  •  physical  quantities,  since  the  coefficients  of 
VQ  and  VQ  in  Equations  (4.3.6)  and  (4.3.7)  have  essentially  no 
analytical  basis.  In  fact,  the  Newtonian  deformation  rates  calculated 
via  Equation  (4.3.9)  (see  Appendix  D)  are  much  larger  than  the  maxi¬ 
mum  possible  deformation  rate  defined  by  Equations  (2.2.4).  Hence, 
in  order  to  preserve  further  development  from  becoming  an  arithmetic 
exercise,  it  is  necessary  to  place  some  severe  restrictions  on  the 
magnitudes  of  the  representative  quantities  . 

As  noted  in  Section  (2.2.1),  the  deformation  rates  of  fibre 
flows  are  expected  to  be  somewhat  lower  than  those  of  Newtonian  flows. 
However,  since  the  actual  deformation  rates  are  unknown,  the  most 
conservative  estimates  will  be  the  "representative"  Newtonian  values 
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given  by  Equation  (4=3.9).  Qualitative  examination  of  the  uncon¬ 
strained  least-squares  coefficients,  [k-j]  and  [k£]  of  all  the 
reliable  tests  performed  in  this  study  (see  Table  (D-l ) )  indicate  the 
following: 

i)  [k-j]  Is  a  decreasing  function  of  the  fibre  concentration 
and  axis  ratio  as  surmised  in  Section  (2.2.2). 
ii)  Numerically,  the  magnitude  of  [k£]  is  very  sensitive  to  the 
magnitude  of  [k-j]. 

Hence,  in  order  to  scale  [k£]  to  the  order  of  the  Newtonian  representa¬ 
tion,  [ k ^ ] ,  a  constrained  least-squares  regression  was  performed  by 


using  the  Newtonian  [k  ]  to  fit  the  fibre  data  in  the  form 


(4.3.10) 


In  this  case,  a  Type  II  fit,  which  weights  points  at  high  velocity, 
best  represented  the  data  and  it  was  used  rather  than  a  Type  III  (see 
Appendix  C) .  The  "representative"  elongational  viscosity.  In),  of 
the  suspensions  can  be  calculated  from 


(4.3.11) 


{n}  =  “D^T  (3w)newt. 


Results  of  these  calculations  are  presented,  in  dejtail,  in  Appendix 
(D)  and  summarized  in  Figure  (4.14). 

The  results  depicted  in  Figure  (4.14)  indicate  that  the 
"representative"  elongational  viscosities  are  much  lower  than  the 
limiting  case  of  Takserman-Krozer  et  al  for  a  =  °°,  which  is  expected 


- 


i.i 


63 


FIGURE  (4.13) 

COMPARISON  OF  "REPRESENTATIVE"  AND 
THEORETICAL  ELONGATIONAL  VISCOSITIES 
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to  apply  for  the  suspensions  used  in  this  study  (see  Section  (1.2.3)). 
However,  the  effects  of  fibre  concentration  and  axis  ratio  is  quali¬ 
tatively  in  accord  with  the  predictions  of  Takserman-Krozer  et  al. 

It  may  also  be  noted  that  the  deformation  rate  used  in  calculating  { n } 
from  the  experimental  results  is  expected  to  be  too  high.  A  less  con¬ 
servative  estimate  may  well  define  elongational  viscosities  which  could 
be  of  the  order  of  those  predicted  by  Takserman-Krozer  and  Ziabicki. 

4.4  UNUSUAL  BEHAVIOUR 

The  unusual  behaviour  mentioned  in  Sections  (4.1)  and  (4.2)  refers 
to  the  apparent  existence  of  two  steady-state  thrust  profiles  with  two 
corresponding  steady-state  pressure  profiles  for  orifice  flows.  Here, 
steady-state  refers  to  a  recorder  trace  of  thrust  or  pressure  which 
remains  constant  with  respect  to  time.  This  effect  appeared  at  random 
in  the  series  of  tests  6  to  16  and  could  not,  in  general,  be  duplicated. 
The  "twin  steady-state"  effect  was  observed  from  a  typical  recorder 
trace  of  pressure  and  thrust  at  constant  flowrate.  This  effect  was 
noted,  as  depicted  in  Figure  (4.14),  for  the  Newtonian  case  in  test  6, 
for  the  polymer  in  test  7,  and  for  the  fibre  suspensions  in  tests  14 
and  16.  An  example  of  this  effect,  when  presented  on  the  basis  of  a 
normalized  thrust,  is  depicted  in  Figure  (4.4)  for  test  16.  The  cases 
where  the  "twin  steady-state"  effect  was  experimentally  encountered 
are  denoted  in  the  results  (see  Appendix  B)  by  two  consecutive  entries 
labelled  "LO"  and  "HI". 
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In  terms  of  measured  thrust,  the  lower  steady-state  value  at  a 
jet  velocity  of  1242  cm./sec.  is  approximately  14%  lower  than  the 
higher  steady-state  value  at  the  same  velocity.  The  difference  in 
0"N)  between  test  9  and  test  16  can  be  explained  on  this  basis.  In 
test  9,  where  the  "twin  steady-state"  behaviour  did  not  appear, 
apparently  only  the  high  values  were  recorded,  whereas,  in  test  16, 
both  effects  were  recorded.  It  is  suspected  that  this  behaviour  may 
have  occurred  in  all  the  tests  from  6  to  16,  but  did  not  display  this 
behaviour  directly  by  providing  two  concurrent  steady-state  pressure 
and  thrust  measurements  during  a  particular  run.  The  results  depicted 
in  Figure  (4.8)  qualitatively  support  this  suspicion.  Occasions 
where  the  unusual  behaviour  was  observed  are  indicated  in  the  data 
(see  Appendix  B)  by  LO  and  HI  denoting  the  low  and  high  values 
respecti vely . 

The  origin  of  the  "twin  steady-state"  effect  for  orifice  flows 
is  presently  not  understood,  hence  the  observations  are  noted.  Since 
this  effect  was  observed  in  tests  6  to  16  where  the  larger  diameter 
(0.1477  cm.)  orifice  was  used,  a  diameter  dependence  is  evident. 

When  the  "twin  steady-state"  effect  occurred,  it  always  went  from 
low  to  high.  Also,  it  may  be  apparent  that  if  a  constant  pressure 
feed  device  is  used  and  flowrates  are  determined  by  sampling  over  a 
time  interval,  an  observation  of  the  "twin  steady-state"  effect  is 
very  unlikely.  These  observations  may  be  useful  in  guiding  future 
experimental  study  of  this  phenomena. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

Thrust  measurements  on  orifice  jets  of  viscous  Newtonian  showed 
abnormal  thrust  reduction  at  low  orifice  velocities.  At  higher 
velocities,  the  viscous  thrust  per  unit  area  asymptotically  approached 
the  thrust  per  unit  area  of  low  viscosity  Newtonians.  An  attempt  to 
relate  the  abnormal  thrust  reduction  to  the  high  tensile  stresses 
developed  in  orifice  flows  was  not  entirely  successful.  Although  the 
profile  relaxation-surface  tension  effect  common  to  capillary  flows 
may  apply  to  orifice  flows,  the  thrust  reduction  for  Newtonians  was 
not  investigated  in  sufficient  detail  to  verify  this  argument. 

The  viscous  Newtonian  thrust  data  were  used  as  a  basis  for 
comparison  with  the  suspension  thrusts.  Orifice  jets  of  fibre  sus¬ 
pensions  exhibited  significantly  lower  thrusts  than  the  viscous 
Newtonian  at  a  given  velocity.  This  thrust  reduction  became  more 
pronounced  with  increasing  axis  ratio  and/or  fibre  concentration. 

If  thrust  reduction  is  regarded  as  an  increase  in  the  axial  stresses 
developed  in  orifice  flows,  the  increased  axial  stress,  due  to  the 
presence  of  fibres,  can  be  related  to  the  elongational  deformation 
rate  by  means  of  an  elongational  viscosity.  Conservative  estimates 
of  the  elongational  viscosities  are  qualitatively  consistent  with  the 
specific  elongational  viscosity  theory  of  Takserman-Krozer  and 
Ziabicki.  Stronger  statements  than  this  are  presently  impossible 
owing  to  the  uncertainty  in  the  velocity  field. 
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Qualitative  observations  have  been  made  of  a  twin  steady-state 
behaviour  in  orifice  flows .  At  present,  no  explanation  is  available 
for  this  phenomena.  Although  the  twin  steady-state  behaviour,  as 
well  as  the  need  for  local  velocity  measurements  complicate  the 
analysis  of  orifice  jet  thrust  measurements,  the  orifice  jet  thrust 
technique  may  become  a  useful  tool  in  non-shearing  rheological 
studies . 

5.2  RECOMMENDATIONS 

The  most  important  recommendations  are  those  directed  towards 
improving  the  equipment  used  in  this  study.  The  suggested  improve¬ 
ments  include  the  following: 

i)  Temperature  sensitivity  of  the  spring  system  could  be 
reduced  by  minimizing  the  bulk  of  the  orifice  reservoir. 

If  this  modification  does  not  reduce  the  temperature 
sensitivity  appreciably,  it  may  be  necessary  to  rebuild 
the  entire  spring  system.  In  case  reconstruction  is 
required,  it  is  suggested  that  only  one  material,  pre¬ 
ferably  mild  steel,  be  used  throughout, 
ii)  Overall  noise  in  the  spring  system  (vibrations)  could  be 
reduced  by  separating  the  spring  leaves  to  a  distance  in 
the  order  of  9  inches.  Small  torsional  displacements  in 
the  axis  of  the  spring  leaves  were  noted  with  the  present 
(heavy)  reservoir  and  leaf  separation  of  5-3/4  inches. 
Addition  of  a  damping  system,  which  is  more  sophisticated 
than  the  present  viscous  dashpot,  would  complete  the  task 
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of  reducing  the  noise  level. 

iii)  A  reservoir  with  a  flush-mounted  pressure  transducer  would 
eliminate  the  present  bother  with  bleeding  the  pressure 
transducer  signal  line. 

The  recommendations  for  future  study  of  suspensions  in  orifice 
flows  are: 

i )  Plugging  or  fibre  aggregation  is  directly  related  to  fibre 
concentration,  fibre  axis  ratio,  and  absolute  particle 
size,  and  inversely  related  to  orifice  size,  suspending 
medium  viscosity  and  fluid  velocity.  Hence,  problems  with 
plugging  or  fibre  aggregation,  at  a  given  fibre  axis  ratio 
and  concentration  can  be  minimized  by  observing  the 
following  guidelines: 

(a)  using  fibres  of  minimum  thickness 

(b)  increasing  orifice  size 

(c)  increasing  viscosity  of  suspending  medium 

(d)  increasing  fluid  velocity 

ii)  The  use  of  sugar  solutions  as  a  suspending  medium  for 

fibres  provides  a  convenient  method  for  suspension  prepara¬ 
tion  since  the  fibres  can  be  easily  dispersed  in  water, 
and  viscosity  can  be  adjusted  by  addition  of  sugar  after 
dispersion  is  complete.  At  high  sugar  concentrations, 
however,  the  suspending  medium  viscosity  is  very  sensitive 
to  small  changes  in  sugar  concentration.  It  is  suggested, 
therefore,  that  viscosity  measurements  be  made  to  ensure 
reproducibi 1 i ty . 


. 

. 


, 


70 


iii)  Since  orifice  flows  of  viscous  media  are  prone  to  a  twin 
steady-state  behaviour,  the  use  of  a  positive  displacement 
pump  is  imperative. 
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FLOWRATE  CALIBRATION 

Since  both  pumping  systems  used  in  this  study  were  positive  dis¬ 
placement  and  could  be  set  at  discrete  speeds,  flow  calibration  was 
a  simple  matter  of  relating  flowrate  to  pump  speed.  For  both  systems, 
calibration  was  accomplished  by  direct  weighing  of  a  fluid  sample 
over  a  known  time  interval.  For  the  Instron  pumping  system,  flowrates 
could  also  be  calculated  since  both  piston  speed  and  cylinder  diameter 
were  accurately  known.  Results  for  the  Zenith  metering  pump  system  are 
tabulated  in  Table  (A-l)  and  plotted  in  Figure  (A-l).  Figure  (A-l) 
shows  that  the  deviation  of  individual  flowrates  from  their  linear 
representation  is  typically  less  than  ±\%.  Calculated  and  calibrated 
flowrates  for  the  Instron  pumping  system  are  tabulated  in  Table  (A-2). 
The  calibrated  flowrates  are  consistently  1%  lower  than  the  calculated 
flowrates,  probably  as  a  result  of  inaccuracies  in  the  calibration 
procedure.  The  calculated  flowrates  were  arbitrarily  used  in  deter¬ 
mining  the  velocities  for  tests  8  to  16  where  the  Instron  pumping 
system  was  used. 
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TABLE  A-l 

ZENITH  PUMP  FLOWRATE  CALIBRATION 


Fluid: 

64%  sugar  solution 

Density: 

Pump 

Speed 

Delivery 

Pressure 

(kd.yne/cm2) 

Fluid 

Mass 

(pi..) 

Time 
(sec. ) 

11.38 

107 

200 

140.3 

14.33 

161 

194 

102.8 

18.04 

226 

198 

87.1 

22.71 

339 

208 

73.1 

28.60 

484 

218 

61.3 

36.00 

709 

204 

45.4 

45.32 

1060 

609 

107.4 

45.32 

1060 

210 

37.0 

57.06 

1610 

658 

92.2 

71  .83 

2380 

699 

77.4 

90.43 

3490 

749 

66.2 

1 .30  gm./cc. 

Flowrate 
(cc./sec. ) 

1.097 

1 .452 

1.749 

2.190 

2.737 

3.458 

4.364 

4.368 

5.492 

6.950 

8.707 


' 

FLOWRATE  (cc/sec) 
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FIGURE  ( A- 1) 

ZENITH  GEAR  PUMP  FLOWRATE  CALIBRATION 


Run 

1 

2 

3 

4 

5 

6 

7 

8 

9 
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TABLE  (A-2) 

INSTRON  FLOWRATE  CALIBRATION 


Fluid:  Water 
Density:  0.997  gm./cc. 


Crosshead 
Speed 
(cm. /min. ) 

Del i very 
Pressure^ 
(kdyne/cmf) 

Fluid 

Mass 

(gm.) 

Time 
(sec. ) 

Measured 

Flowrate 

(cc./sec) 

Calculated 

Flowrate 

(cc/sec) 

10 

0 

3374 

112.0 

30.21 

30.40 

5 

0 

3033 

201.2 

15.12 

15.20 

2 

0 

1108 

184.6 

6.02 

6.08 

1 

0 

1124 

374.6 

3.01 

3.04 

7 

0 

3167 

150.1 

21.16 

21.28 

3.5 

0 

2082 

198.0 

10.54 

10.64 

1.4 

0 

1247 

296.7 

4.24 

4.26 

7 

7000 

2515 

119.0 

21.19 

21  .28 

5 

7000 

2500 

165.4 

15.16 

15.20 
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PRESSURE  TRANSDUCER  CALIBRATION 


Pressure  transducer  calibration  was  accomplished  by  comparing 
amplifier  output  of  the  pressure  transducer  signal,  to  manometer 
reading.  In  order  to  cover  the  entire  range  of  pressures  encountered 
experimentally,  two  manometers  were  used;  one  containing  mercury 
(S.G.  =  13.53)  (30),  and  one  containing  water  (S.G.  =  0.997)  (30). 
Tables  (A-3)  and  (A-4)  represent  the  mercury  and  water  data  respect¬ 
ively.  These  results,  which  represent  three  ranges  in  transducer 
response,  appear  graphically  in  Figures  (A-3),  (A-4)  and  (A-5).  For 
each  range,  the  pressure  can  be  expressed  as  a  linear  function  of 
amplifier  output.  These  linear  expressions  and  their  ranges  of 
applicability  are  as  follows: 


P  =  5.215  V 


P  =  5.05  V 


P  =  9.00  +  4.96  V 


(0  <  V  <  40) 

(40  <  V  <  180) 

(180  <  V  <  600)  (A-l ) 


2 

where  P  is  the  pressure  in  kdyne/cm  .  and  V  is  amplifier  output  in 


arbitrary  units. 


' 

. 


Run 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 
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TABLE  A-3 

PRESSURE  TRANSDUCER  CALIBRATION 


Pressure  Voltage 

(inches  Hg)  (kdyne/cm  . )  (arbitrary  units) 


79.0 

2665 

528 

76.8 

2590 

512 

72.8 

2460 

486 

67.2 

2270 

454 

61.9 

2090 

414 

59.1 

1993 

394 

57.0 

1923 

382 

48.0 

1620 

324 

44.6 

1503 

300 

41.5 

1400 

280 

39.6 

1337 

268 

37.4 

1262 

252 

35.4 

1194 

238 

33.4 

1127 

223 

29.4 

992 

196 

27.2 

918 

183 

24.5 

826 

165 

23.35 

787 

157 

22.30 

753 

150 

19.60 

661 

133 

18.10 

611 

119 

16.80 

567 

112 

15.65 

528 

103 

12.85 

434 

85.5 

12.40 

418 

82.5 

11.30 

381 

75.5 

10.35 

349 

68.5 

9.25 

312 

62.0 

8.10 

273 

52.8 

7.65 

258 

49.6 

6.30 

212 

41.4 

5.50 

185.4 

36.6 

4.60 

155.2 

30.4 

3.95 

133.2 

25.6 

2.65 

89.5 

17.2 

' 

• 

Run 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


TABLE  A-4 
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PRESSURE  TRANSDUCER  CALIBRATION 


Pressure  2  Voltage 

(inches  H^O )  (kdyne/cm.  j  (arbitrary  units) 


3.20 

7.96 

1.51 

6.10 

15.18 

2.84 

10.25 

25.5 

4.72 

13.4 

33.3 

6.35 

15.4 

38.4 

7.30 

22.4 

55.6 

10.50 

26.2 

65.0 

12.2 

31.3 

77.9 

14.9 

38.3 

95.4 

18.3 

50.4 

125.3 

23.6 

58.4 

145.2 

28.6 

69.4 

173 

33.4 

73.3 

182 

35.0 

78.4 

195 

37.4 

80.5 

200 

38.6 

PRESSURE  (kdy  ne/cm2) 
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0  10  20  30  40  50  60 


VP- VOLTAGE  (arb.  units) 

FIGURE  (A- 3) 

PRESSURE  TRANSDUCER  CALIBRATION 


PRESSURE  (kdyne/cm2) 
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FIGURE  (A-4) 

PRESSURE  TRANSDUCER  CALIBRATION 


PRESSURE  (kdyne/cm?) 
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FIGURE  (A- 5) 

PRESSURE  TRANSDUCER  CALIBRATION 
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THRUST  CALIBRATION 

Thrust  calibration  was  accomplished  by  relating  amplifier  out¬ 
put  of  the  displacement  transducer  signal  to  a  known  force  applied 
normal  to  the  spring  face  and  in  the  axis  of  the  orifice.  The  known 
force  was  applied  through  the  use  of  a  suspended  balance  weight,  WT, 
as  depicted  schematically  in  Figure  (A-6).  For  a  given  WT,  the  hori¬ 
zontal  force,  FH,  is  given  by 


FH  =  (0.981)  (WVHWTj  (A_2) 

-  (WV)2 

where  WT  is  expressed  in  gm.,  FH  in  kdyne,  and  WV  and  H  in  cm.  The 
thrust  calibrations  were  not  reproducible  from  test  to  test  due  to 
the  sensitivity  of  the  spring  system  to  small  variations  in  mounting 
the  reservoir.  Hence,  a  separate  calibration  was  required  for  each 
test  in  this  study. 

An  estimate  of  the  maximum  possible  error  in  calculating  the 
applied  force,  FH,  is  given  by 


(AFH)  •  £ 


3FH 

3WT 


jAWTj 


3FH 

3WV 


AWV 


aFH 


3H 


ah 


(A-3) 


where  the  quantities  prefixed  by  A  indicate  the  error  in  that  quantity 
and 

3 FH  (0.981)  (WV) 

3WT  '  rH2-(WV)¥;"Z 

3FH  =  (0.981)(MT)(H2) 

3WV  [H2-(WV)2]3/2 


I 


' 

'to  s$£t  tP9  nA 


i 
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FIGURE  (A-6) 

THRUST  CALIBRATION  SCHEME 
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3FH  =  -(0,981 HWTHWVHhO 

3H  [H2-(WV}2]3/2 


(A-4) 


Arbitrarily  choosing  run  2  of  the  thrust  calibration  for  test  9  (see 
Table  A-13),  as  a  typical  measurement, 


WT  =  5.452  gm. 

W V  =  21.77cm. 

H  =  39.20  cm. 

FH  =  3.572  kdyne 

Assuming  that  the  error  in  WT  is  the  thread  mass,  which  was  found  to 
be  approximately  0.06  gm.,  and  the  errors  in  W V  and  H  to  be  the  small¬ 
est  division  in  the  rule  (0.1  cm.),  the  error  in  FH  will  be 

aFH  =  (0.655) (0.06)  +  (0 .237) (0 . 10)  +  (0 . 1 31 3) (0 . 1 ) 

=  0.0761  kdyne 

A  typical  (maximum)  error,  therefore,  is  0.0761  kdyne  or  2.1%  of  the 
applied  force.  This  error  is  smaller  than  the  error  in  reading  the 
transducer  signal  on  the  recorder,  which,  among  other  things,  is 
affected  by  the  zero  drift  of  the  spring  system. 

The  use  of  two  different  orifice  sizes  in  this  study  required 
a  change  in  the  overall  sensitivity  of  the  transducer-carrier  preamp 
system  (or  a  change  in  the  Cal  Factor)  to  accommodate  the  maximum 
sensitivity  of  the  transducer  for  each  orifice  size.  Hence,  in 
order  to  express  the  orifice  thrust  in  terms  of  an  amplifier  output 
which  would  be  equivalent  for  both  sensitivities  (or  Cal  Factors), 
an  intermediate  variable  was  defined,  namely: 


■ 


' 
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VSTAR  =  Amplifier  output  in  volts/Cal  Factor  ( A-5 ) 

The  thrust,  T,  could  then  be  represented  as, 

T  =  B  VSTAR  (A-6) 

where  the  constant  B  was  determined  from  a  first  order  TYPE-II 
least-squares  regression  of  FH  versus  VSTAR. 

Results  of  the  thrust  calibrations  are  presented  in  Tables 
(A-5)  to  (A-21)  along  with  their  respective  least-squares  representa¬ 
tions.  An  example  of  a  calibration  change  during  a  run  is  presented 
in  Tables  (A-18)  and  (A-19)  which  are  a  before  and  after  calibration 
for  test  14.  Although  the  difference  in  the  least-squares  representa¬ 
tions  is  small  (less  than  1%),  the  individual  measurements  of  thrust 
are  expected  to  vary  from  ±3%  at  the  higher  values  to  ±10%  at  the 
lower  values  of  measured  thrust.  This  error  is  expected  to  result 
mainly  from  the  temperature  sensitivity  of  the  spring  system. 


- 
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TABLE  (A-5) 


THRUST  CALIBRATION  FOR  TEST  NO.  1 


LENGTH  OF  HYPOTENUSE  =  32.930  CM. 
CALIBRATION  FACTOR  *  1.000 
VSTAR  =  VT/CAL  FACTOR 


WT 

WV 

VT 

FH 

VSTAR 

(GM) 

(CM) 

(VOLTS  AT  XI) 

( KDYNES ) 

4.452 

17.920 

26.400 

2.833 

26.400 

5.452 

17.550 

32.100 

3.369 

32.100 

4.452 

22.950 

42.000 

4.245 

42.000 

5.452 

22.010 

48.600 

4.806 

48.600 

7.452 

20.850 

63.100 

5.981 

63.100 

9.452 

19.990 

70.800 

7.084 

70.800 

9.452 

26.760 

122.300 

12.931 

122.300 

14.452 

24.910 

162.200 

16.399 

162.200 

THRUST (KDYNE )  =  0.10183  *  VSTAR 
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TABLE  (A-6) 


THRUST  CALIBRATION  FOR  TEST  NO.  2 


LENGTH  OF  HYPOTENUSE  =  32.930  CM. 
CALIBRATION  FACTOR  *  1.030 
VSTAR  =  VT/CAL  FACTOR 


WT 

(GM) 

WV 
(  CM) 

VT 

(VOLTS  AT 

FH 

XI)  (KDYNES) 

VSTAR 

4.452 

17.380 

30.400 

2.714 

29.514 

5.452 

17.090 

36.200 

3.247 

35.145 

4.452 

22.760 

47.400 

4.177 

46.019 

5.452 

21.960 

53.000 

4.787 

51.456 

7.452 

20.430 

68.500 

5.733 

66.504 

9.452 

19.850 

82.000 

7.006 

79.611 

14.452 

19.540 

108.500 

10.453 

105.339 

10.452 

26.400 

158.000 

13.754 

153.398 

14.452 

24.770 

186.000 

16.186 

180.582 

19.452 

23.360 

222.000 

19.209 

215.533 

24.452 

22.490 

256.000 

22.431 

248 . 543 

THRUST (KDYNE ) 

=  0.09023 

*  VSTAR 

' 
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TABLE  (A-7) 

THRUST  CALIBRATION  FOR  TEST  NO.  3 


LENGTH  OF  HYPOTENUSE  =  33.000  CM. 
CALIBRATION  FACTOR  *  1.030 
VSTAR  =  VT/CAL  FACTOR 


WT 
(  GM) 

WV 

(CM) 

VT 

(VOLTS  AT 

FH 

XI)  ( KDYNES ) 

VSTAR 

4.452 

18.030 

24.600 

2.849 

23.683 

5.452 

17.620 

31.600 

3.378 

30.679 

4.452 

22.940 

41.200 

4.224 

40.000 

6.452 

21.250 

50.500 

5.328 

49.029 

9.452 

19.800 

67.000 

6.955 

65.048 

14.452 

18.640 

94.000 

9.706 

91.262 

9.452 

26.760 

125.000 

12.851 

121.359 

14.452 

24.820 

161.000 

16.182 

156.310 

24.452 

22.410 

214.000 

22.194 

207.766 

29.452 

21.560 

243.000 

24.937 

235.922 

THRUST (KDYNE ) 

=  0.10588 

*  VSTAR 

1  :  J 
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TABLE  (A-8) 

THRUST  CALIBRATION  FOR  TEST  NO.  4 

LENGTH  OF  HYPOTENUSE  =  33.000  CM. 
CALIBRATION  FACTOR  *  1.030 
VSTAR  =  VT/CAL  FACTOR 

WT  WV  VT  FH  VSTAR 

(GM)  (CM)  (VOLTS  AT  XI)  ( KDYNES ) 


4.452 

17.530 

25.200 

2.738 

24.466 

5.452 

17.150 

32.600 

3.253 

31.650 

4.452 

22.910 

43.000 

4.213 

41.747 

7.452 

20.810 

65.500 

5.940 

63.592 

9.452 

19.940 

72.000 

7.032 

69.9C2 

14.452 

18.480 

96.500 

9.584 

93.689 

19.452 

17.900 

122.500 

12.322 

118.932 

14.452 

24.890 

168.000 

16.288 

163.106 

24.452 

22.320 

221.000 

22.030 

214.563 

29.452 

21.550 

249.000 

24.917 

241.747 

THRUST (KDYNE ) 

=  0.10215  * 

VSTAR 

. 
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TABLE  (A-9) 

THRUST  CALIBRATION  FOR  TEST  NO*  5 

LENGTH  OF  HYPOTENUSE  =  33.000  CM. 
CALIBRATION  FACTOR  *  1.030 
VSTAR  =  VT/CAL  FACTOR 


WT 

(GM) 

WV 

(CM) 

VT 

(VOLTS  AT  XI) 

FH 

( KDYNES ) 

VSTAR 

4.452 

22.810 

40.000 

4.178 

38.634 

5.452 

21.870 

45.200 

4.734 

43.683 

7.452 

20.760 

57.500 

5.917 

55.825 

9.452 

19.870 

67.000 

6.993 

65.048 

14.452 

18.560 

93.000 

9.645 

90.291 

19.452 

17.730 

125.500 

12.157 

121.844 

14.452 

24.750 

158.000 

16.078 

153.398 

24.452 

22.480 

214.000 

22.324 

207.766 

29.452 

21.700 

241.000 

25.222 

233.980 

THRUST (KDYNE ) 


0.10638  *  VSTAR 
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TABLE  (A-10) 

THRUST  CALIBRATION  FOR  TEST  NO.  6 

LENGTH  OF  HYPOTENUSE  =  33.000  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  *  VT/CAL  FACTOR 


W  T 
(GM) 

WV 

(CM) 

VT 

(VOLTS  AT  XI) 

FH 

( KDYNES ) 

VSTAR 

4.452 

22.740 

21.200 

4.153 

39.186 

5.452 

21.930 

24.000 

4.757 

44.362 

9.452 

19.760 

35.400 

6.933 

65.434 

14.452 

18.590 

46.600 

9.667 

86.136 

14.452 

24.820 

84.000 

16.182 

155.268 

19.452 

23.340 

98.500 

19.094 

182.070 

24.452 

22.310 

112. 000 

22.012 

207.024 

29.452 

21.650 

131.000 

25.119 

242.144 

34.452 

20.980 

145.000 

27.841 

268.022 

54.452 

19.480 

202.000 

39.071 

373.382 

THRUST (KDYNE  ) 


0.10472  *  VSTAR 
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TABLE  (A-l 1 ) 


THRUST  CALIBRATION  FOR  TEST  NO.  7 

LENGTH  OF  HYPOTENUSE  =  33.000  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

WV 

VT 

FH 

VSTAR 

( GIY ) 

(CM) 

(VOLTS  AT 

XI)  ( KDYNES ) 

4.452 

22.740 

19.700 

4.153 

36.414 

5.452 

21.930 

22.700 

4.757 

41.959 

7.452 

20.780 

28.600 

5.926 

52.865 

9.452 

19.760 

33.600 

6.933 

62.107 

14.452 

18.590 

45.600 

9.667 

84.288 

14.452 

24.820 

78.500 

16.182 

145.101 

24.452 

22.310 

106.000 

22.012 

195.933 

34.452 

20.980 

136.000 

27.841 

251.386 

54.452 

19.480 

187.000 

39.071 

345.656 

THRUST (KDYNE ) 

=  0.11227 

*  VSTAR 

. •••  - 
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TABLE  (A- 12) 


THRUST  CALIBRATION  FOR  TEST  NO.  8 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

(GM) 

WV 

(CM) 

VT 

(VOLTS  AT  XI) 

FH 

( KDYNES ) 

VSTAR 

4.452 

23.380 

17.800 

3.245 

32.9C2 

5.452 

23.050 

21.200 

3.888 

39.186 

6.452 

22.840 

24.600 

4.538 

45.471 

9.452 

22.360 

35.400 

6.440 

65.434 

5.452 

28.250 

29.600 

5.560 

54.713 

9.452 

25.980 

44 .400 

8.207 

82.070 

14.452 

24.550 

62.500 

11.390 

115.526 

19.452 

23.720 

79.000 

14.505 

146.025 

24.452 

23.200 

96.000 

17.615 

177.449 

24.452 

28.690 

141.000 

25.767 

260.628 

34.452 

27.200 

177.000 

32.572 

327.171 

54.452 

25.430 

249.000 

45.542 

460.258 

104.452 

23.540 

426.000 

76.965 

787.430 

114.452 

23.330 

460.000 

83.165 

850.277 

THRUST (KDYNE ) 

=  0.09813  *  VSTAR 

• 

. 

• 
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TABLE  (A-l 3) 


THRUST  CALIBRATION  FOR  TEST  NO.  9 


LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

(GM) 

WV 

(CM) 

VT 

(VOLTS  AT  XI) 

FH 

( KDYNES ) 

VSTAR 

4.452 

23.880 

22.900 

3.355 

42.329 

5.452 

21.770 

28.400 

3.572 

52.495 

5.452 

28.210 

38.600 

5.544 

71.349 

9.452 

25.970 

57.500 

8.202 

106.284 

14.452 

24.540 

79.000 

11.383 

146.025 

24.452 

23.200 

119.500 

17.615 

220.887 

24.452 

28.690 

178.000 

25.767 

329.020 

34.452 

27.180 

223.000 

32.526 

412.199 

54.452 

25.420 

320.000 

45.511 

591.497 

104.452 

23.520 

534.000 

76.863 

987.060 

THRUST (KDYNE ) 


0.07784  *  VSTAR 


►  >’ 
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TABLE  (A-14) 

THRUST  CALIBRATION  FOR  TEST  NO.  10 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


W  T 

WV 

VT 

FH 

VSTAR 

(  GM  ) 

(CM) 

(VOLTS  AT  XI) 

( KDYNES ) 

4.452 

23.330 

21.600 

3.235 

39.926 

5.452 

22.950 

25.800 

3.663 

47.669 

5.452 

28.210 

37.200 

5 . 544 

68.761 

9.452 

25.960 

55.500 

8.196 

102.587 

14.452 

24.520 

76.000 

11.368 

140.480 

24.452 

23.190 

119.500 

17.603 

220.887 

34.452 

22.530 

162.000 

23.741 

299.445 

24.452 

28.710 

177.000 

25.806 

327.171 

34.452 

27.220 

222.000 

32.618 

410.351 

54.452 

25.440 

310.000 

45.573 

573.012 

104.452 

23.540 

522.000 

76.965 

964.879 

94.452 

23.790 

480.000 

70.763 

887.245 

THRUST (KDYNE  )  =  0.07967  *  VSTAR 
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TABLE  (A-15) 


THRUST  CALIBRATION  FOR  TEST  NO.  11 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

WV 

VT 

FH 

VSTAR 

(GM) 

(  CM) 

(VOLTS  AT  XI) 

( KDYNES ) 

4.452 

23.390 

22.700 

3.247 

41.959 

5.452 

22.980 

26.800 

3.870 

49.537 

5.452 

28.220 

37.400 

5.548 

69.131 

9.452 

25.960 

59.000 

8.196 

109.057 

14.452 

24.530 

81.000 

11.375 

149.722 

24.452 

23.200 

123.000 

17.615 

227.356 

34.452 

22.560 

168.000 

23.788 

310.536 

24.452 

28.650 

182.000 

25.690 

336.414 

34.452 

27.130 

230.000 

32.411 

425.138 

54.452 

25.400 

328.000 

45.449 

606.284 

84.452 

24.050 

460.000 

64.376 

850.277 

104.452 

23.530 

548.000 

76.914 

1012.939 

THRUST (KDYNE ) 


0.07583  #  VSTAR 


• 

' 

•  »C  I 

* 
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TABLE  (A- 16) 

THRUST  CALIBRATION  FOR  TEST  NO.  12 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

(GM) 

WV 

(CM) 

VT 

(VOLTS  AT  XI) 

FH 

( KDYNES ) 

VSTAR 

4.452 

23.360 

23.100 

3.241 

42.698 

5.452 

22.940 

27.600 

3.860 

51.016 

5.452 

28.260 

39.000 

5.564 

72.088 

9.452 

25.980 

59.000 

8.207 

109.057 

14.452 

24.550 

81.000 

11.390 

149.722 

24.452 

23.190 

123.500 

17.603 

228.280 

34.452 

22.550 

171.000 

23.772 

316.081 

24.452 

28.720 

183.000 

25.826 

338.262 

34.452 

27.210 

230.000 

32.595 

425.138 

54.452 

25.450 

330.000 

45.604 

609.981 

84.452 

24.080 

484.000 

64.505 

894.639 

104.452 

23.550 

552.000 

77.016 

1020.332 

THRUST (KDYNE )  =  0.07452  *  VSTAR 
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TABLE  (A-17) 


THRUST  CALIBRATION  FOR  TEST  NO.  13 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

WV 

VT 

FH 

VSTAR 

(GM) 

(CM) 

(VOLTS  AT 

XI)  ( KDYNES ) 

4.452 

23.360 

19.200 

3.241 

35.489 

5.452 

22.950 

22.900 

3.863 

42.329 

5.452 

28.410 

33.800 

5.626 

62.476 

9.452 

25.960 

48.600 

8.196 

89.833 

14.452 

24.300 

69.000 

11.201 

127.541 

24.452 

23.190 

105.500 

17.603 

195.009 

34.452 

22.560 

145.000 

23.788 

268.022 

24.452 

28.700 

167.000 

25.787 

308.687 

34.452 

27.190 

197.000 

32.549 

364.140 

54.452 

25.440 

280.000 

45.573 

517.560 

84.452 

24.050 

394.000 

64.376 

728.280 

104.452 

23.540 

464.000 

76.965 

857.670 

THRUST (KDYNE ) 

=  0.08875 

*  VSTAR 

. 
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•TABLE  (A-18) 

THRUST  CALIBRATION  FOR  TEST  NO.  14 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

(GM) 

WV 

(CM) 

VT 

(VOLTS  AT  XI) 

FH 

( KDYNES ) 

VSTAR 

4.452 

23.390 

20.200 

3  •  247 

37.338 

5.452 

22.960 

23.700 

3.865 

43.807 

5.452 

28.190 

35.000 

5.536 

64.695 

9.452 

25.940 

53.500 

8.185 

98.890 

14.452 

24.520 

72.500 

11.368 

134.011 

24.452 

23.200 

109.500 

17.615 

202.402 

34.452 

22.560 

150.000 

23.788 

277.264 

24.452 

28.680 

163.000 

25.748 

301.293 

34.452 

27.160 

205.000 

32.480 

378.927 

54.452 

25.520 

294.000 

45.821 

543.438 

84.452 

24.050 

410.000 

64.376 

757.855 

104.452 

23.520 

484.000 

76.863 

894.639 

THRUST(KDYNE)  *  0.08537  *  VSTAR 
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TABLE  (A-19) 


THRUST  CALIBRATION  FOR  TEST  NO.  14 

LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

WV 

VT 

FH 

VSTAR 

(GM) 

(CM) 

(VOLTS  AT  XI) 

( KDYNES ) 

4.452 

23.330 

20.600 

3.235 

38.077 

5.452 

22.940 

23.300 

3.860 

43.068 

5.452 

28.260 

34.600 

5.564 

63.955 

9.452 

25.980 

49.800 

8.207 

92.051 

14.452 

24.540 

71.500 

11.383 

132.162 

24.452 

23.180 

110.500 

17.591 

204.251 

34*452 

22.540 

152.000 

23.756 

280.961 

24.452 

28.700 

165.000 

25.787 

304.990 

34.452 

27.200 

207.000 

32.572 

382.624 

54.452 

25.440 

296.000 

45.573 

547.134 

84.452 

24.090 

408.000 

64.548 

754.158 

104.452 

23.540 

486.000 

76.965 

898.336 

THRUST (KDYNE )  =  0.08520  *  VSTAR 


• 

. 
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TABLE  (A-20) 

THRUST  CALIBRATION  FOR  TEST  NO.  15 


LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

WV 

VT 

FH 

VSTAR 

(GM) 

(CM) 

(VOLTS  AT  XI) 

( KDYNES ) 

4*452 

23.350 

21.000 

3.239 

38.817 

5.462 

22.980 

24.300 

3.870 

44.916 

5.452 

28.220 

36.200 

5.548 

66.913 

9.452 

25.970 

54.500 

8.202 

100.739 

14.452 

24.540 

76.000 

11.383 

140.480 

24.452 

23.200 

114.000 

17.615 

210.720 

34.452 

22.540 

156.000 

23.756 

288.354 

24.452 

2e.700 

171.000 

25.787 

316.081 

34.452 

27.200 

213.000 

32.572 

393.715 

54.452 

25.410 

306.000 

45.480 

565.619 

84.452 

24.070 

428.000 

64.462 

791.127 

104.452 

23.520 

508.000 

76.863 

939.001 

THRUST (KDYNE ) 


0.08164  #  VSTAR 


■ 

p'i 
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TABLE  (A-21) 

*  \  ■ 

THRUST  CALIBRATION  FOR  TEST  NO.  16 


LENGTH  OF  HYPOTENUSE  =  39.200  CM. 
CALIBRATION  FACTOR  *  0.541 
VSTAR  =  VT/CAL  FACTOR 


WT 

wv 

VT 

FH 

VSTAR 

(GM) 

(CM) 

(VOLTS  AT 

XI)  ( KDYNES ) 

4.452 

23.290 

20.900 

3.226 

38.632 

5*452 

22.900 

24.000 

3.850 

44.362 

5.452 

28.160 

35.600 

5.523 

65.804 

9.452 

25.940 

51.400 

8.185 

95.009 

14.452 

24.500 

73.000 

11.352 

134.935 

24.452 

23.150 

112.000 

17.556 

207.024 

34.452 

22.510 

155.000 

23.709 

286.506 

24.452 

28.640 

168.000 

25.671 

310.536 

34.452 

27.120 

211.000 

32.368 

390.018 

54.452 

25.380 

300.000 

45.386 

554.528 

84.452 

24.050 

418.000 

64.376 

772.643 

104.452 

23.490 

498.000 

76.710 

920.517 

THRUST (KDYNE ) 

=  0.08309 

*  VSTAR 

' 

" 

. 
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APPENDIX  B 

EXPERIMENTAL  RESULTS 


B-l 


THRUST  MEASUREMENTS 

Experimental  results  are  summarized  in  Tables  (B-l)  to  (B-16). 

The  heading  of  each  table  classifies  the  fluid  according  to  the  weight 
concentration  of  sugar,  axis  ratio  and  weight  concentration  of  fibres, 
and  gives  the  orifice  size  used.  Occasions  where  an  unusual  behaviour 
(see  Section  4.4)  occurred  are  indicated  by  two  consecutive  entries 
at  the  same  velocity  marked  LO  and  HI  to  the  right  of  the  column 
headed  by  "TN".  LO  and  HI  denote  the  low  and  high  values,  respectively, 
which  characterize  the  "unusual  behaviour". 
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APPENDIX  C 

SAMPLE  CALCULATIONS  AND 
LEAST-SQUARES  REGRESSIONS 
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SAMPLE  CALCULATIONS 


All  calculations,  including  the  least-squares  regressions, 
were  performed  with  the  aid  of  the  IBM  1800  Real-Time  computer  of  the 
DACS  Centre.  The  sample  calculations  presented  below  are  based  on 
the  first  run  of  test  9. 

Recorder  signals  of  thrust  and  pressure  were  converted  to  the 
physical  units  given  in  Tables  (B-l)  to  (B-16)  by  the  procedure  out¬ 
lined  in  Appendix  A.  From  Table  (B-9) 


Q  =  3.040  cc./sec. 


T  =  0.627  kdyne 
P  =  97.52  kdyne/cm.^ 


The  velocity  of  the  orifice,  VQ,  is  given  by 


V 


3.040  cc/sec 


177.428  cm. /sec. 


o 


The  thrust  per  unit  orifice  area  is  given  by 


T  T 


0.627  kdyne  =  36-616  |<dyne/cm2 

,0.1477/ 


The  normalized  thrust,  TN,  defined  by  Equation  (4.2.2),  is 


TN 


36.616 


=  1.858 


(1 ,8667)(10'3)(177.428)2  -  (0 . 22014) (177.428) 


The  averaged  normalized  thrust,  (TN),  defined  by  Equation  (4.2.3),  is 


■ 
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(1 .68)(10~3)(200)2  -  ( 0 . 1 89 ) ( 200 ) 

(1 .8667)(10"3)(200)2  -  (0.22014) (200) 


0.959 


for  an  orifice  velocity  of  200  cm. /sec.  The  contraction  coefficient, 
C  ,  is  defined  by  Whitaker  (44),  and  for  a  large  contraction  ratio 

V/ 

(ratio  of  reservoir  area  available  for  flow  to  orifice  area)  is 
approximated  by 


Cc  -  ^ 

For  the  first  run  of  test  9,  Cc  is 

Cc  =  177.428  / (2) (97ilf (.1000)  =  0,461 


The  results  of  the  above  calculations  are  given  in  Tables  (B-l)  to 
(B-16). 


( 
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LEAST-SQUARES  REGRESSION 

If  a  function  y  s  f(x)  can  be  adequately  represented  by  a  poly¬ 
nomial,  then  the  well-known  principle  of  least-squares  (38)  can 
easily  be  applied  to  obtain  the  "best-fit"  polynomial  which  defines 
the  function  y.  For  example,  data,  (x^ ,  y^ ) ,  which  can  be  expressed 
as  a  second-order  polynomial  in  x,  will  give  the  following  system  of 
"residual"  equations: 


r 

L 

2 

n  EX..  EX.j 

ao 

2  3 

EX.  EX.  EX.'3 

al 

= 

Exiyi 

2  3  4 

2 

EXi  EX.j  EXi 

_a2_ 

[_Exi  yiJ 

or,  in  matrix  notation, 

X  A  =  Y  (C-2) 

where  x^,(i  =  l,2,3,...,n)  are  the  data  points  of  the  independent 
variable  x,  and  y^ ,  (i  =  1,2,3,...,n)  are  the  corresponding  data 
points  of  the  dependent  variable,  y.  The  solution  to  the  system  of 
Equations  (C-l)  is  given  by 

y  =  a0  +  a-|X  +  a2x2  (C-3) 

such  that  the  sum  of  the  squares  of  the  deviations  of  y^,(i  =  1,2,3, 
. ,„,n)  from  y  will  be  a  minimum.  For  convenience,  the  unconstrained 
least-squares  polynomial  defined  by  Equation  (C-3)  will  be  labelled 
as  a  Type  I  fit. 


.... 


. 
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If  apriori  knowledge  of  the  function  y  indicates  that  y(0)  =  0, 
then  Equation  (C-2)  can  be  constrained  to  force  the  coefficient  aQ  in 
Equation  (C-3)  to  zero  without  violating  the  least-squares  principle. 
At  least  two  approaches  are  possible: 

1.  Apriori  specification  of  aQ  =  0  leads  to  the  following 
system  of  residual  equations  (38) 


2 

EX.j 

EX.j 

al 

EVi 

3 

4 

ii 

2 

LExi 

ZX.j 

r\3 

1 _ 

_Exi  yi_ 

(C-4) 


The  least-squares  fit  derived  from  Equation  (C-4)  will  be  labelled  as 
a  Type  II  fit, 

2o  A  simple  constraint  is  provided  by  transforming  Equation 
(C-3)  to 


z  =  x  =  ai  +  a2x  (C-5) 

giving  a  Type  I,  first-order  polynomial  whose  solution 
can  be  obtained  from  the  "residual"  system  defined  as 


n 

M 

X 

_ 1 

af 

N 

_ J 

EX- 

M 

X 

ro 

1 

_a2_ 

EX.Z. 

_  l 

! 

_  1  1_ 

(C— 6 ) 


The  least-squares  fit  available  from  the  solution  of  (C-6)  will  be 
labelled  as  a  Type  III  fit. 

A  solution  to  either  system  of  Equations  (C-4)  or  (C-6)  will 
result  in  a  polynomial  of  the  form 


. 


. 
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y  =  a-jX  +  a2x2  (C-7) 

which  requires  that  y  pass  through  the  origin.  The  difference  between 
a  Type  II  and  a  Type  III  fit  can  be  qualitatively  described  on  the 
basis  of  the  difference  in  the  order  of  x  in  the  corresponding 
elements  of  (C-4)  and  (06).  Since  the  order  of  x  in  an  element  of 
(04)  is  +2  greater  than  the  order  of  x  in  a  corresponding  element 
of  (06),  a  Type  II  fit  will  be  more  sensitive  to  higher  values  of  x 
than  a  Type  III  fit.  Basically,  a  Type  II  fit  will  provide  a  better 
fit  at  essentially  high  values  of  x,  while  a  Type  III  fit  may  provide 
a  better  overall  fit. 

A  statistical  goodness-of-fi t  criterion,  such  as  the  correla¬ 
tion  coefficient  (4),  will  generally  indicate  that  a  Type  I  fit 
represents  data  better  than  a  Type  II  or  a  Type  III  fit  since  the 
Type  I  fit  is  less  constrained  than  the  other  two.  However,  apriori 
knowledge  of  a  functional  form  often  suggests  a  constraint  which  is 
based  (from  a  statistical  point  of  view)  on  personal  judgement  and 
is  not  taken  into  account  by  the  goodness-of-fi t  criterion.  Hence, 
goodness-of-fi t  was  determined  by  personal  judgement  rather  than 
statistical  means. 

A  qualitative  comparison  of  the  Type  II  and  Type  III  fits  is 
given  in  Figures  (C-l)  and  (C-2).  Figure  (C-l)  depicts  the 
"unconstrained"  thrust-velocity  data  representation  of  test  13  by 
each  fit.  Figure  (C-2)  depicts  the  "constrained"  data  representation 
of  test  13  by  each  fit.  It  is  apparent  that  a  Type  III  fit  best 
represents  the  "unconstrained"  thrust-velocity  relationship  while 


L  I 


•  r  V  ill  .  .  , 

. 
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a  type  II  fit  best  represents  the  "constrained"  version. 
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FI6URE  (C-1) 

UNCONSTRAINED  REPRESENTATION  OF  TEST  13 


-  . 
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FIGURE  (C-2) 

CONSTRAINED  REPRESENTATION  OF  TEST  13 
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ELONGATIONAL  VISCOSITIES 


1  ■  * 


D-l 


"REPRESENTATIVE"  ELONGATIONAL  VISCOSITY 


A  more  or  less  qualitative  prediction  of  the  elongational 
viscosity  behaviour  of  fibre  suspensions  can  be  made  via  the  develop¬ 
ment  presented  in  Section  (4.3.2).  Accordingly,  the  least-squares 
representation  of  the  thrust-flowrate  data  for  the  sugar  solutions 
were  used,  in  conjunction  with  the  least-squares  coefficients  given 
in  Table  (D-l)  to  predict  the  "representative"  elongational  viscosity, 
{n>,  of  each  suspension  used  in  the  reliable  tests  of  this  study. 

Since  the  suspending  medium  for  test  3  was  a  60%  sugar  solu¬ 
tion,  [k^]  is  given  by  Equation  (4.2.4),  namely: 

[k2]  =  0.1639 

The  shear  viscosity,  y,  of  the  60%  sugar  solution  is  (43) 


u  =  0.52  poise 

The  fibre  concentration,  given  in  Table  (3.1)  and  the  "constrained" 

least-squares  coefficient,  [k,!>R],  given  in  Table  (D-l)  are 

C  =  0.090  vol .  % 

[4r]  =  0.221 


Hence,  for  test  3,  In)  is  given  by  Equation  (4.3.11),  i.e. 


{n}  =  ( 3) (0 . 52 )  =  2.10  poise 

In  order  to  find  what  the  elongational  viscosity,  { n >66-2/3%  wou^c*  be 
for  a  volumetric  concentration  of  0.0926%,  a  fibre  axis  ratio  of 


. 


85  and  suspended  in  66-2/3%  sugar  solution,  Equation  (1.2.6)  is 
utilized,  i .e. 


D-2 


2.10  -  1.56 _ 

(9 .00) (10~4) (1 .56) 


385 


Since 


n 


then 


{n}66-2/3%  =  5,1[1  +  (0.000926) (385)] 


=  6.92  poise 

It  is  now  possible  to  compare  -C n ^66-2/3%  e^on9atlona^ 


viscosities  of  the  suspensions  used  in  tests  9,  11,  13,  and  16. 

The  "representative"  elongational  viscosities  of  the  sus¬ 
pensions  used  in  the  reliable  tests  of  this  study  were  calculated 
as  outlined  above  and  appear  in  Table  (D-2).  Table  (D-2)  also 
lists  the  theoretically  predicted  (39,40)  elongational  viscosities, 

n '  and  n ' . 

o  °° 


THEORETICAL  ELONGATIONAL  VISCOSITY 


Rudd  (34)  defines  a  rotational  diffusion  constant,  D,  which 
is  given  by 


3kT[1  n(^jj-)  -  1] 

75 

*  ^o1 


(D-l) 


D 


where  k  is  the  Boltzman  constant,  T  the  absolute  temperature, 


« 


:•  ’  « 


D-3 


the  shear  viscosity  of  the  suspending  medium,  and  1  and  d  are  the 

length  and  diameter,  respectively,  of  the  particle.  For  the  fibres 

-12  -1 

used  in  this  study,  D  has  a  value  of  10  sec.  Since 


a 


(D-2) 


the  limiting  case  of  a  =  °°  is  expected  to  apply  for  the  suspensions 
used  in  this  study. 

Equation  (13)  of  Takserman-Krozer  and  Ziabicki  (40)  define 
v  for  the  limiting  cases  of  a  =  0  and  a  =  °°  for  a  suspension  of  rigid 
rods ,  i .e. 


where 

K 

M 

N 


=  ^(2K  +  7M  +  12N) 
=  (2K  +  M) 


1 


1 

b;[(ar)2  +  i] 


(D-3) 


(D-4) 


~a'Q  =  |o[2(AR)2  -  5  +  fx] 


For  rods, 


« 


D-4 


ii  — 

a  -  aLo 

o  u2 


ll  (AR)2  +  l  +  iLAg|-— 1  x] 


s:  =  +i  - 1  x] 

0  (AR)2  2 


6 


“  =  0[-3  ♦  Mliil  x] 


(D-5) 


where 


a  = 


1 _ ln  AR  f  /(AR)2  -  1 


AR/(AR)2  -  1  AR  -  /(AR)2  -  1 


(0-6) 


The  intrinsic  viscosities,  \>  and  v  ,  were  calculated,  as  outlined 

o  00 

above,  for  the  fibres  used  in  this  study  and  the  "limiting"  elonga- 

tional  viscosities,  n1  and  n1  were  determined  from  Equation  (1.2.9), 

0  00 

namely: 


u  =  (“C^C  -  0  (D'7) 

and  Equation  (1 .2.6) 

-  n'  ~  °0 

nsp  n0  (0-8) 

The  limiting  elongational  viscosities  and  n^,  for  the  suspensions 
used  in  tests  9,  10,  11,  13,  14  and  16  are  presented  in  Table  (D-2). 


' 

. 


■ 


■ 


D-5 


TABLE  D-1 

LEAST-SQUARES  COEFFICIENTS  OF 

RELIABLE  THRUST-VELOCITY  DATA 


Test 

[k]] 

3 

0.00163 

9 

0.00168 

10 

0.00160 

11 

0.00161 

13 

0.00156 

14 

0.00163 

16 

0.00162 

[k^]  [k£R] 


0.100 

0.221 

0.189 

0.361 

0.162 

0.418 

0.191 

0.433 

0.221 

0.515 

0.130 

0.363 

0.164 

0.402 

D-6 


TABLE  D-2 

"REPRESENTATIVE11  AND  THEORETICAL 

ELONGATIONAL  VISCOSITIES 


Test 

AR 

Fibre 

Cone. 

(n> 

n'o 

i 

— 

(wt.  %) 

(poise) 

(poise) 

(poise) 

3 

85 

0.10 

6.92 

6.95 

14.7 

9 

170 

0.10 

8.35 

11.40 

36.6 

10 

170 

0.20 

9.65 

16.1 

78.4 

11 

255 

0.10 

10.00 

18.2 

70.0 

13 

340 

0.10 

11.90 

26.9 

137.0 

14 

85 

0.20 

8.37 

9.08 

23.9 

16 

170 

0.10 

9.27 

11.4 

36.6 

